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ABSTRACT 

An  attempt  to  describe  the  atmospheric  boundary  layer  during 
the  MIZEX  84  from  11  July  to  17  July  was  made  using  a  sodar 
system.  The  computation  of  the  temperature  structure  parameter  CT^ 
in  the  surface  layer  using  in  situ  measurements  allowed  the  calibra- 
tion of  the  system.  CT^  was  found  to  have  the  following  functional 
dependence  on  the  backscattered  signal  I: 
Log(CT2)=    8.63Log(I)     -25.66. 

Time-height  cross-sections  showed  very  well  the  structure  of 
the  boundary  layer.  An  especially  good  representation  of  the 
strength    of    CT^    in    the    inversion    layer    was    achieved. 

A  program  was  developed  to  obtain  the  inversion  layer  height 
Zi,  where  the  signal  reaches  a  maximum.  The  thickness  of  the 
inversion  layer  was  also  computed  and  was  used  with  the  calibration 
law  to  compute  CT^  at  the  inversion  layer.  This  allowed  the 
computation  of  the  jump  of  potential  temperature  and  the  refractive 
index  gradient  in  the  inversion  under  free  convection.  Refractive 
trapping  conditions  were  evident  on  15  July  although  the  generally 
small  value  of  the  heat  flux  did  not  indicate  a  strongly  unstable 
boundary  layer.  Comparaison  with  two  radiosonde  profiles  shows 
good  agreement  in  one  case,  but  does  not  permit  a  reliable  conclu- 
sion to  be  made  concerning  the  method.  However,  the  procedure 
could  be  extended  to  other  areas,  especially  those  with  strong  free 
convection.  The  method  should  provide  a  good  estimate  of  the 
refractive    condition    in    the    atmospheric    boundary    layer. 


RESUME 

Un  systeme  sodar  (Sound  Detection  and  Ranging)  a  permis  de 
recueillir  des  donnees  concernant  la  couche  limite  atmospherique 
entre  le  11  et  17  juiilet  pendant  I'operation  MIZEX  1984.  Le  calcul 
du  parametre  de  structure  de  la  temperature  CT^  dans  la  couche  de 
surface  utilisant  des  mesures  in  situ,  permet  la  calibration  du 
systeme.  II  apparait  que  CT^  est  fonction  du  signal  de  reverbera- 
tion   direct    I    suivant    la    loi: 

logCCT^)    =    log(I)     -    25.66. 

Plusieurs  coupe  representant  CT^  en  fonction  du  temps  et  de 
I'altitude  montrent  la  structure  de  la  couche  limite  et  donnent  une 
bonne  representation  de  la  valeur  de  CT^  dans  la  couche 
d'inversion. 

Un  programme  a  ete  developpe  pour  obtenir  la  hauteur  de  cette 
couche  o.u  un  signal  .maximum  est  atteint.  L'epaisseur  de  cette 
couche  est  aussi  calculee  afin  d'etre  utilisee  avec  la  loi  de  calibra- 
tion pour  en  deduire  CT^  a  chaque  niveau.  Ceci  permet  de 
trouver  le  saut  de  temperature  potentieUe  et  le  gradient  d'index  de 
refraction  dans  I'inversion  en  cas  de  convection  libre.  Pendant  la 
journee  du  15  juUIet  on  remarque  que  les  ondes  electromagnetiques 
ont  pu  etre  "capturees"  par  la  couche  d'inversion,  bien  que  la 
faible  valeur  du  flux  de  chaleur  n'ait  pas  ete  toujours  conforme  au.x 
conditions      d'instabilite      requises.  La      comparaison      avec      quelques 

profiles  de  radiosonde  est  bonne  dans  un  cas,  mais  ne  permet  pas 
d'obtenir  une  conclusion  certaine  sur  la  methode.  II  pourrait  cepen- 
dant  etre  interessant  de  I'etendre  a  d'autres  zones  soumises  a  une 
plus  forte  convection  libre,  ceci  donnant  une  estimation  des  condi- 
tions   de    refraction    dans    la    couche    limite    atmospherique. 
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I.    INTRODUCTION. 


The  atmospheric  boundary  layer  (ABL),  which  is  found  just 
above  the  surface  and  generally  extends  to  several  hundred  meters, 
can  be  considered  to  be  composed  of  three  different  layers.  At  the 
lowest  level  is  the  surface  layer,  whose  height  is  about  one  tenth 
of  the  ABL  height  and  where  turbulent  exhanges  with  the  surface 
are  driven  directly  by  vertical  gradients.  The  mixed  layer  comprises 
most  of  the  remaining  ABL  depth.  Under  unstable  conditions 
mechanical  and  thermal  mixing  cause  homogenous  conditions.  The 
third  layer  is  the  inversion  layer  where  turbulent  kinetic  energy  is 
small  and  there  are  strong  vertical  gradients.  In  particular,  the 
gradient  of  potential  temperature  is  associated  with  stronger  fluctua- 
tions   of    the    temperature. 

The  polar  region  is  of  particular  interest  because  it  is  a 
region  not  only  of  scientific  interest  but  also  of  military  importance. 
This  area  differs  significantly  from  the  midlatitude  and  tropical 
ABL's,  where  anomalous  electromagnetic  propagation  due  to  mean 
vertical  gradient  of  temperature  and  humidity  are  frequently 
observed.  Although  anomalous  propagation  is  not  commonly  reported 
in  the  Arctic  Ocean,  the  boundary  layer  can  extend  to  a  low 
height  and  have  a  strong  inversion  layer.  Another  difference 
concerns  the  difference  in  the  stability  condition,  with  a  generally 
more  stable  atmosphere  at  high  latitudes  and  more  convective  condi- 
tions in  mid-latitude  or  in  the  tropics.  Therefore,  it  is  important 
to  understand  the  vertical  structure  of  Arctic  atmospheric  boundary 
layer.        To    date,    little    work    have    has    been    done    on    this    problem. 

The  Marginal  Ice  Zone  experiment  (MIZEX  84)  was  conducted  in 
the  summer  of  1984  and  produced  a  wealth  of  data  with  which  to 
investigate  the  structure  of  the  Arctic  ABL.  Several  instrumented 
platforms     were      deployed     in     the      East     Greenland     Sea      during     that 
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summer  to  measure  the  mesoscale  processes  that  influence  the 
advance  and  retreat  of  the  ice  margin.  A  variety  of  data  were  then 
gathered,  including  measurements  of  surface  variables  and  sodar 
data.  The  sodar  (Sound  Detection  and  Ranging)  system  is  one  of 
the  sensors  used  in  this  study.  It  measures  the  strength  of  tiie 
small-scale  turbulence  encountered  along  its  vertical  ray  path.  The 
intensity  of  the  signal  is  proportional  to  the  refractive-index  struc- 
ture parameter  Cn^,  called  also  the  structure  function  parameter, 
which  has  an  important  environmental  influence  on  the  propagation 
of  electromagnetic  waves.  Profiles  of  this  parameter  provided  by  the 
sodar  allow  production  of  an  image  of  the  structure  of  the  atmos- 
pheric   boundary    layer. 

It  is  the  purpose  of  this  thesis  to  calibrate  the  sodar  system 
deployed  by  the  research  vessel  POLAR  QUEEN  during  MIZEX  84  to 
obtain  the  refractive -index  structure  parameter.  The  time-height 
structure  of  the  boundary  layer  will  also  be  investigated.  Efforts 
will  also  be  made  to  compute  the  height  of  the  inversion  layer  and 
to  estimate  its  thickness,  the  jump  of  potential  temperature,  and  the 
refractive-index    gradient. 
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II.    THE   EXPERIMENT 

A .       LOCATION 

1.  Overview 

The  recent  phase  of  the  Marginal  Ice  Zone  Experiment  was 
conducted  during  the  summer  of  1984.  This  experiment,  which 
followed  that  of  1983,  was  designed  to  explore  the  drift,  the  forma- 
tion, and  the  evolution  of  the  ice  margin  through  mesoscale 
processes. 

The  experiment  began  with  a  preparation  period  consisting 
mainly  of  satellite  imagery  analysis.  Following  this  period,  seven 
ships,  eight  aircraft  and  four  helicopters  belonging  to  the  ten 
european  and  north -american  countries  involved  in  the  experiment 
were  deployed  in  Fram  Strait  area  between  Greenland  and  Svalbard 
Island  (Spitsbergen).  The  first  platform  was  on  station  at  the 
beginning  of  June  and  the  recovery  period  started  on  17  July 
(Johannessen    and    Horn,     1984) . 

2.  Polar   Queen    Station 

The  main  goal  of  the  POLAR  QUEEN  was  to  foUow  the 
drifting  motion  of  the  ice.  Therefore,  she  was  moored  to  a  floe  a 
few  kilometers  inside  the  ice  edge  and  drifted  passively  with  the 
ice  floe  while  various  experiments  were  deployed  to  study  the 
parameters  which  described  the  atmospheric  boundary  layer.  Two 
different  drifting  phases  occurred.  The  first  one  ended  on  16  June 
when  the  floe  was  broken  by  the  propagating  swell  generated  by 
northerly  winds.  The  station  was  then  moved  about  60  kilometers  to 
the    northwest    and    observations    terminated    on    17    July. 

This  most  important  part  of  the  experiment  was  located  in 
an    area     delineated    by    latitudes    80°N     to    81°N    and     longitudes    001°E 
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to  007°E.  We  will  limit  this  study  to  the  last  week  of  the  experi- 
ment from  11  July  at  2150Z  to  17  July  at  1130Z  where  a  most 
interesting  data  set  was  obtained.  During  this  period  the  sodar 
display      terminal      showed       a      well      formed      inversion      layer.  This 

terminal  produced  a  time-height  cross  section  of  the  backscattered 
signal  using  shaded  characters  whose  strength  is  proportional  to  the 
signal  intensity.  This  gave  a  good  in  situ  picture  of  the  .ABL 
structure. 


Fdllback   Floe,    1   July   1984,   Scale   1:7200 


Figure    2.1  The    Polar    Queen    Station. 

The  ship  was  moored  on  the  port  side  to  a  large  floe  of 
approximately  400  by  700  meters  as  shown  in  figure  2.1.  The 
POLAR    QUEEN    can    be     seen    in    dark,     moored    along    the    upper    edge 
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of  the  main  floe  (named  "Fallback  floe").  As  it  can  be  seen  in 
this  picture,  the  surrounding  region  is  composed  of  floes  of  a 
variety  of  sizes.  The  floes  are  all  drifting,  but  not  necessarily  with 
the  same  motion.  So  the  environment  of  the  station  varied,  and  this 
can  induce  variations  in  the  heat  and  moisture  flux  depending  on 
location.  A  parcel  of  air  which  is  at  a  low  height  will  be  under 
the  influence  of  the  underlying  ocean  surface.  A  higher  parcel  can 
be    affected    by    both    ice    and    water    surfaces. 


B.       WEATHER    CONDITIONS 

The  weather  was  characterized  generally  by  stratus  or  fog 
conditions,  although  there  was  less  fog  than  during  the  1983  exper- 
iment. Several  low  pressure  systems  passed  through  the  region 
during  the  experiment.  For  the  period  of  interest,  the  area  was 
characterized  by  a  high  pressure  system  to  the  east  of  Svalbard 
Island.  Another  weaker  high  pressure  center  lay  over  Greenland, 
while  a  low  pressure  system  was  located  north  of  Greenland  to  the 
northwest  of  the  POLAR  QUEEN  station.  These  systems  moved 
slowly  and  erratically.  The  pressure  varied  from  1004  millibars  to 
1020  millibars,  with  the  minimum  values  occurring  on  about  11  July, 
with  higher  values  later.  No  front  crossed  the  area  during  that 
week.  Winds  were  variable  with  velocities  ranging  from  0  to  15 
m/sec. 


C .       MEASUREMENTS 

The    two    different    sources    of    the    data    to    deal    with    were    the    in 
situ    measurements    and    the    sodar    measurements. 
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1.  In    Situ    Measurements 

In  situ  measurements  were  obtained  at  two  different  loca- 
tions. One  was  the  ship  itself,  where  surface  pressure,  air 
temperature,  dew  point  temperature,  wind  speed  and  relative  direc- 
tion,   and    downward    radiation    were    measured. 

The  second .  location  was  an  ice  mast  located  on  the  ice 
floe  at  about  50  meters  from  the  edge  as  shown  on  figure  2.1  at 
spot  "A".  It  provided  air  temperature  and  wind  speed  at  different 
levels,     from    5    cm    to    5.3    m    above    the    floe    surface. 

2.  Sodar   Measurements 


To  avoid  noise  interference  from  the  ship,  the  SODAR 
system  was  installed  on  the  floe  at  site  "B"  as  is  shown  in  figure 
2.1,  about  30  meters  from  the  ice  edge.  It  is  expected  that  the 
acquired  measurements  were  more  influenced  by  surface  ice  than  by 
the    surrounding    water. 

The  sodar  system  will  be  described  in  more  detail  in  the 
next  chapter.  However,  the  system  gave  essentially  a  backscattered 
signal  which  is  proportional  to  the  small  scale  fluctuation  of  temper- 
ature in  a  layer  encountered  by  the  acoustic  signal.  These  fluctua- 
tions influence  the  refractive  index  structure  parameter  Cn^,  which 
is  mainly  dependent  upon  the  temperature  structure  parameter  CT^. 
This  parameter  is  generally  larger  in  the  surface  layer  and  in  the 
inversion  (or  interfacial)  layer.  The  sodar  gave  the  values  of  tlie 
backscatterred  signal  for  each  height  at  a  rate  of  about  one  profUe 
every  minute.  It  will  be  possible  to  deduce  the  height  of  the 
inversion  layer,  and  some  information  regarding  its  strength  may  be 
inferred.  One  goal  of  this  study  is  to  calibrate  the  sodar  by 
computing  the  value  of  CT^  in  the  surface  layer  with  only  in  situ 
measurements  and  then  comparing  these  values  with  the  sodar  signal 
strength    to    obtain    a    calibration    curve. 
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D.        GENERAL   DATA    REDUCTION. 

1.       Determination       of       Surface      Layer       Scales       Obtained      from 
Surface    Layer  ~^imiTarity 

The  in  situ  measurements  include  humidity  and  wind 
speed  and  temperature  at  different  levels.  From  these  U=",  T'"  and 
q"'  may  be  obtained  using:  U^==<-u'w'>^  ^^,  T"'=<-w'T'>/U^'=  and 
q=!-=<-w'q'>/U'''-.  The  bracket  represents  an  average  over  time. 
These  quantities  were  computed  using  the  bulk  method  (Businger, 
1973)      derived      from      Monin-Obukhov      similarity      theory.  A      brief 

description    of    the    method    and    results    follows. 

According  to  the  bulk  method,  the  following  equations  can 
be    used    to    compute    U'''    and    T-": 

U"''    =    kU<z>    /    (ln(Z/Zo)     -    0,(O)  (2.1) 

T='=    =    ka(T-To)    /    (ln(Z/Zo)     -    0    (C))  (2.2) 

2 

where: 

k    =    0.35  von    Karman's    constant. 

a    =    1.35  ratio    of    heat    transfer    to    momentum    transfer    at    ^    =    0. 

U<z'     (m/s)       wind    speed    at    height    Z. 

Zo,     Zot  roughness    length    for    velocity    and    temperature    profiles. 

0     (C)    and     0     (^)     closed    functions    of    the     similarity    height    parameter 

^  is  defined  as  the  ratio  Z/L  where  L  is  the 
Monin-Obukhov    length    scale    (Monin    and    Obukhov,     1954) 

E,    =    Z/L    =    Zkg    (T-    +6.1    lO-^Tci'-')    /    T    U^:^^  (2.3) 

where     g    =     9.8    m/s^     is    the    acceleration    of     gravity    and    q'''     (gm/kg) 
is    the    water    vapor    mixing    ratio    scaling    parameter. 

5  was  computed  neglecting  q'^'  because  of  the  high  relative 
humidity  and  low  temperature  which  produce  a  low  humidity  gradient 
and,     hence,     a    low    value    for    q=''.     Equation    2.3    becomes: 
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C    =    Zkg    T'^    /    T    U=i='  (2.4) 

Equation  2.1  can  be  put  in  the  form  U<z'  =  m  f(Z,0  *  ^,  where 
f(Z,C)  can  be  computed  using  an  iteration  process  to  get  C-  m  = 
U'Vk  is  the  slope  of  the  curve  and  can  be  measured  giving  the 
value  of  U'''.  n  =  U'''ln(Zo)/k  is  the  intercept  and  allows  the  compu- 
tation of  Zo.  In  fact,  the  drag  coefficient  was  first  computed  using 
the    formula : 

Cd^/2    =    U=VU<'''>  (2.5) 

where  U*^°'  is  the  wind  speed  at  10  meters.  The  range  of  "Cd" 
was    found    to    be    between    2.1    and    2.3    10"^    m/s. 

T='=    is    compute    in    the    same    way    as    U''=    and    the    two    results 
produce    the    following    transfer    coefficient: 

Ch    =    U=!=T-    /    U<z'(T-To) 

where  T=V(T-To)  =  0.032  is  obtained  from  a  plot  of  T='=  versus 
(T-To)    This    leads    to: 

Ch    =    1.5    10"^ 

(Personal    communication    with    P.     Guest,     1985). 


iRg 


vironmental    Considerations    for    the    Calculation    of    the    Heat 


The  heat  flux  can  be  computed  easily  once  we  know  U""" 
and    T- 

Qo    =    -    U-T-  .  (2.G) 

A  summary  of  the  most  useful  data  for  this  study  obtained 
from  the  in  situ  measurements  and  averaged  over  15  values  (about 
every    3    hours)     are    presented    in    Table    I  Here    the    character     "'-^" 

means    that    there    are    no    reliable    data    for    the    corresponding    period. 
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TABLE   I 
IN   SITU   AND    REMOTE   MEASUREMENTS 


Date  -  Time 

Temp 

Rh. 

m/s  . 

P 
mbar. 

U"       T" 
m/s       °K 

Qo 
°K(m/s) 

llJul 

0453Z 

0.3 

99.1 

4.8 

1004.0 

0.2270 

0.0166 

-0.0038 

llJul 

0730Z 

0.5 

99.8 

5.6 

1004. 1 

0.2622 

0.0205 

-0.0054 

llJul 

1002Z 

0.3 

98.8 

5.0 

1004.5 

0.2229 

0.0120 

-0.0027 

llJul 

1236Z 

0.2 

99.0 

5.9 

1004.9 

0.2697 

0.0094 

-0.0025 

llJul 

1509Z 

0.0 

98.1 

5.5 

1005.3 

0.2504 

0.0046 

-0.0012 

llJul 

1741Z 

-0.2 

97.3 

4.7 

1005.7 

0.2065 

-0.0013 

0.0003 

llJul 

2027Z 

-0.6 

95.2 

3.4 

1006.7 

0.1055 

-0.0086 

0.0009 

llJul 

2259Z 

-0.4 

90.2 

2.7 

1007.9 

0.0653 

0.0029 

-0.0002 

12Jul 

0131Z 

-0.6 

89.6 

2.1 

1009.0 

0.0583 

-0.0002 

0.0000 

12Jul 

0402Z 

-0.8 

93.0 

1.2 

1010.0 

0.0359 

-0.0389 

0.0014 

12Jul 

0634Z 

-0.9 

89.6 

1.4 

1011.2 

0.0618 

-0.0411 

0.0025 

12Jul 

0906Z 

-1.2 

89.3 

1.8 

1012.2 

0.0529 

-0.0864 

0.00h6 

12Jul 

1137Z 

-1.0 

89.5 

2.3 

1012.9 

0.0626 

-0.1853 

0.0116 

12Jul 

1411Z 

-1.3 

92.3 

2.8 

1013.3 

0.0384 

-0.1944 

0.0075 

12Jul 

1654Z 

-1.5 

94.8 

3.2 

1013.4 

0.0210 

-0.0766 

0.0016 

12Jul 

1927Z 

-1.2 

96.6 

3.8 

1013.3 

0. 1165 

-0.0434 

0.0051 

12Jul 

2159Z 

-0.7 

99.7 

3.8 

1013.2 

0. 1009 

-0.0164 

0.0017 

13Jul 

0032Z 

-0.4 

100.0 

3.3 

1012.7 

0.1577 

-0.0124 

0.0019 

13Jul 

0304Z 

-0.1 

100.6 

2.2 

1012.1 

0.0309 

-0.0250 

0.0008 

13Jul 

0535Z 

-0.1 

100.7 

2.4 

1011.5 

0.0488 

-0.0175 

0.0009 

13Jul 

0811Z 

-0.1 

100.6 

2.8 

1011.3 

0.0765 

-0.0205 

0.0016 

13Jul 

1042Z 

-0.3 

99.6 

5.2 

1011.4 

0.1181 

-0.0432 

0.0051 

13Jul 

1313Z 

-0.2 

99.3 

6.2 

1011.2 

0.2978 

-0.0303 

0.0090 

13Jul 

1544Z 

-0.5 

99.  1 

5.6 

1011.4 

0.2681 

-0.0327 

0.0088 

13Jul 

1814Z 

0.5 

100.2 

3.9 

1011.5 

9.:J19A 

-0.0060 

0.0016 

13Jul 

2045Z 

1.0 

100.5 

4.8 

1011.8 

0.0047 

-'.",'.--.'.--;.- v.- v.--;,- 

13Jul 

2317Z 

1.2 

99.9 

5.0 

1011.6 

0,'0689 

0.0017 

-q.oppj. 

14Jul 

0250Z 

1.1 

99.4 

4.2 

1011.4 

*'*»*»  * r  *' v*" 

0.0032 

14Jul 

0522Z 

0.4 

99.4 

4.4 

1011.5 

_0,2358 

-0.0024 

,o,pop6 

14Jul 

0754Z 

0.8 

99.4 

4.8 

1011.7 

0.0023 

14Jul 

1025Z 

0.5 

98.1 

5.3 

1011.7 

2.0735 

-0.0319 

0.0023 

14Jul 

1257Z 

0.6 

98.5 

5.3 

1011.7 

-0.0082 

"Jc  ~  V  Vc  V'  Vr  Vr  *  V 

14Jul 

1529Z 

0.6 

99.8 

4.1 

1011.6 

0.0219 

-0.0090 

'0.0002 

14Jul 

1803Z 

-0.2 

99.8 

6.0 

1011.4 

0.1264 

-0.0131 

0.0017 

14Jul 

2034Z 

-0.1 

99.9 

6.3 

1011.5 

0. 1368 

-0.0019 

0.0003 

14Jul 

2306Z 

-0.0 

100.1 

6.3 

1011.2 

0.0872 

0.0016 

-0.0001 

15Jul 

0138Z 

-0.2 

99.3 

6.6 

1011.2 

0.2367 

-0.0015 

0.0003 

15Jul 

0409Z 

-0.4 

95.3 

7.2 

1011.2 

0.3067 

-0.0095 

0.0029 

15Jul 

0641Z 

-0.1 

92.9 

7.4 

1011.6 

0.3020 

0.0015 

-0.0005 

15Jul 

0913Z 

0.0 

99.5 

6.8 

1012. 1 

0.2702 

0.0019 

-0.0005 

15Jul 

1144Z 

-0.2 

98.0 

6.6 

1012.4 

0.2972 

-0.0081 

0.0024 

15Jul 

1416Z 

-0.3 

93.7 

7.0 

1012.8 

0.3306 

-0.0117 

0.0039 

15Jul 

1741Z 

-0.6 

97.9 

6.4 

1013.3 

0.3024 

-0.0205 

0.0062 

15Jul 

2102Z 

-0.4 

94.8 

7.0 

1013.8 

0.3140 

-0.0072 

0.0023 

15Jul 

2334Z 

0.1 

99.  1 

6.2 

1014.4 

0.2397 

0.0164 

-0.0039 

16Jul 

0207Z 

0.5 

97.9 

5.9 

1014.7 

0.2544 

0.0281 

-0.0072 

16Jul 

0439Z 

0.8 

98.6 

5.4 

1015.0 

0.2316 

0.0297 

-0.0069 

16Jul 

0746Z 

1.2 

99.8 

4.7 

1015.6 

0.1991 

0.0408 

-0.0081 

16Jul 

1017Z 

1.3 

99.9 

4.7 

1016.3 

0.2077 

0.0420 

-0.0087 

16Jul 

1248Z 

1.5 

98.5 

5.3 

1016.7 

0.2104 

0.0486 

-0.0102 

16Jul 

1519Z 

1.5 

99.3 

5.4 

1016.9 

0.1910 

0.0486 

-0.0093 

16Jul 

1750Z 

1.7 

99.2 

6.4 

1017.0 

0.1804 

0.0487 

-0.0088 

16Jul 

2021Z 

1.7 

99.4 

6.8 

1017.3 

0,2056 

.0.0572 

-0,0117 

16Jul 

2251Z 

1.7 

99.8 

7.4 

1017.8 

^CrCfC^C  "i-C^yC 

17Jul 

0122Z 

1.9 

99.5 

7.9 

1018.1 

iVVf-.VVf-;.--;'.- 

-,v -;.--;.- VfVr -A- -iV 

-;.-  -;.-  -;.-  -.v  -;.-  v.-  -.v 

17Jul 

0354Z 

2.2 

98.2 

8.5 

1018.3 

""'''"'"'''''"' 

^C  *C ^C*'*»  *'^C 

•>- .'-  .••  •'-  .•-  -•-  .•- 
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A  question  which  arises  at  this  point  is:  "What  was  the 
relative  influence  of  the  sea  and  of  the  ice  on  the  computation  of 
U"-''  and  T'''  ?"  For  the  ice  mast,  there  was  nearly  no  influence  of 
sea  water  because  the  measurements  are  taken  close  to  the  floe 
surface.  However  for  the  ship,  which  was  in  tlie  water  surrounded 
by  scattered  pieces  of  ice,  the  problem  may  be  more  complex.  The 
surface  of  the  ice  can  de  considered  to  be  nearly  at  the  same 
temperature  as  the  sea  water,  close  to  273°K.  It  can  also  be 
assumed  that  the  humidity  of  the  ice  surface  was  nearly  equal  to 
100  %  because  of  a  generally  thin  skin  of  melting  ice  and  snow  on 
the  floes,  which  more  or  less  resembled  large  ponds  of  water  in 
this  last  part  of  the  experiment.  Therefore,  there  will  be  no 
effort  to  separate  the  influences  from  water  and  from  ice  in  the 
derivation    of    U^^^,     T='=,     q='=    and    Qo. 

3.       Calculation    of    q* 

Assuming  than  the  previous  coefficient  Ch  can  be  applied 
to  the  formula  which  gives  q"%  it  is  possible  to  e.xamine  the 
assumption  that  the  effects  of  q=''  are  negligible.  An  alternative 
way    to    express    q'''    and    T'"    is    the    following: 

rj.,,    3    ^1/2     (x-To)  (2.7) 

CV'-    =    c'^^     (q-qo)  (2.8) 

where  "q"  and  "qo"  are  the  specific  humidity  at  10  meters  and  at 
the  surface  respectively,  "c"  is  the  drag  coefficient  and  is  related 
to    Ch    by: 

c^^^    =    Ch    U<^°VU-  (2.9) 

If  the  drag  coefficients  are  the  same  for  q-  and  T-'-,  q-"'  can  be 
expressed    by: 
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q-    =    1.5    10-'    U<'°>(q    -    qo)/    U*  (2.10) 

(q  -  qo)  is  computed  using  Clausius-Clapeyron  equation  (Fleagle 
and    Businger,     1980) : 

q-qo=    exp[(9.4-2353/T)LnlO](Rh/100    -1)0.622/P  (2.11) 

Here  P  is  the  pressure  in  millibars  and  Rh  is  the  relative 
humidity  on  the  ship.  The  measurement  is  assumed  to  apply  to  a 
height  of  10  meters.  The  relative  humidity  at  the  surface  was  taken 
to  be  100%  because  of  the  large  ponds  of  water  on  the  ice  floe  and 
the    melting    snow    as    stated    previously. 

The  magnitude  for  T-  was  between  10"^  °K  and  10"^  °K 
while  6.1  10"'*Tq='=  goes  from  10"^  to  10"^  g/kg-  So  this  term  can, 
in    fact,     be    neglected    in    equation    2.3    . 

4.       The    Roughness    Length 
Assuming    that    Cd    is    measured    close    to    neutral    stability    we    have: 

Cd'^2    =    k/    ln(Z/Zo)  (2.12) 

which    leads    to:  • 

Zo    =    Z    exp(-k/    Cd'/^'*  (2.13) 

for    Z    =    10    meters    and    Cd    =    2.2    10"',     we    get: 

Zo  =  2  millimeters.  This  value  is  greater  than  the  usual 
quantity  for  the  roughness  length,  probably  because  of  the 
irregularity    of    the    ice    surface. 
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III.    SODAR    OPERATION    AND    ANALYSIS 

A.       DESCRIPTION    OF    THE   SYSTEM 
1.       Basic    Principles 

The  acoustic  sounder  operates  much  like  a  classical  radar 
but  emits  an  acoustic  rather  than  an  electromagnetic  pulse  in  a 
narrow  beam  and  receives  a  signal  which  comes  from  scattering  of 
the  acoustic  wave  by  air  parcels  at  various  levels.  This  backscat- 
tered  signal  is  due  to  small-scale  acoustic  refractive  index  inliomo- 
geneities  which  are  the  consequences  of  local  fluctuations  of  the  air 
temperature  and  velocity.  The  effects  of  the  humidity  fluctuations 
are  negligible  because  they  do  not  affect  significantly  the  sound 
velocity  fluctuations,  which  are  directly  related  to  the  changes  in 
index  of  refraction  (Parry  et  al.  ,  1972).  This  will  be  verified 
later , 

A  monostatic  sodar  system  uses  an  antenna  configuration  'in 
which  the  transmitters  and  receivers  are  co-located.  Thus,  the 
received  signal  is  due  to  direct  backscatter.  Since  the  velocity 
variations  produce  a  maximum  scattering  in  the  forward  direction 
and  do  not  contribute  to  backscattering,  the  received  intensity  will 
not  depend  on  the  air  velocity  fluctuations  and  will  be  mainly  a 
function  of  the  local  temperature  variations.  These  are  described  by 
the  temperature  structure  parameter  "cx^".  Since  temperature  fluc- 
tuations are  largest  in  the  surface  layer  and  in  the  inversion  layer 
of  an  unstable  ABL,  there  will  be  maximum  backscattered  signal 
intensity    in    these    regions. 
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2.       Characteristics    of    the    System 

The  system  emits  a  100  ms  pulse  at  a  frequency  f=  2  kHz 
and  then  shifts  to  reception  mode  for  a  few  seconds.  This  time  is 
variable,  depending  on  the  height  being  probed.  The  listening 
period  is  divided  into  several  segments  or  "bursts"  of  4  ms,  each 
which  will  represent  a  certain  distance  covered  by  the  signal.  To 
have  a  complete  representation  of  the  vertical  layer  above  the 
antenna,  the  total  sounding  depth  is  divided  into  altitude  intervals 
called  "gates"  which  are  characterized  by  a  certain  number  (n)  of 
bursts.  An  average  is  made  for  each  gate  to  give  the  final  backs- 
cattered  signal.  The  altitude  interval  d  of  the  gates  is  manually 
variable  by  changing  the  number  n.  The  value  of  d  is  then  given 
by    the    formula : 

d    =    4    10"'    n    V    /2,  (3.1) 

where  V  is  the  speed  of  sound  in  the  atmosphere  and  is  approxi- 
mated   by    the    formula: 

V    =    20    T'^2    m/s 

where  T  is  the  absolute  temperature  of  the  medium.  For  T  =  273°K 
we  get  V  =  331  m/s,  which  was  used  to  determine  the  height 
scale.  Consequently,  for  a  gate  of  7  bursts,  the  interval  will  be 
d=  4.6  meters.  There  are  200  gates  for  each  profile  which  can 
cover  a  height  of  25  +  (200  4.6)  =  945  meters  which  will  be  a 
maximum  height  of  the  backscattered  profiles.  The  measurements 
begin  at  25  meters  because  of  transmitter  ringing  and  antenna 
near-field    effects. 

Another  parameter  of  interest  is  the  wavelength  of  the 
signal    which    is    given    by: 

L    =    V/f    =    16.5    cm 

The  most  important  scattering  by  small  scale  fluctuations  is  done  by 
eddies  whose  size  is  half  the  wave  length  (about  8  cm  in  our  case) 
of    the    acoustic    wave    (Parry    et    al.  ,     1975). 
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3.       Dependence    of    the    Received    Signal    on    the    Layer    Stability 

The  backscattered  signal  depends  on  the  small  scale  inho- 
mogeneities  encountered,  which  are  related  to  the  atmospheric 
boundary    layer    structure    and    stability. 

a.  Unstable    Mixed    Layer 

Hydrostatic  instability  occurs  when  the  surface  temper- 
ature is  larger  than  the  temperature  of  the  layer  above.  This 
induces  upward  (positive)  buoyancy  flux  (and  heat  flux  Qo) .  Tlie 
resultant  mixing  extends  to  the  interfacial  layer.  Then  the  poten- 
tial temperature  away  from  the  surface  becomes  more  homogeneous 
and  its  variations  with  height  are  reduced,  producing  less  backscat- 
tering.  The  highest  values  of  the  backscatter  signal  level  will  be  in 
the  surface  layer  and  in  the  interfacial  layer  where  the  vertical 
gradients    of    temperature    are    much    greater. 

b.  Stable    Boundary    Layer. 

When  the  surface  temperature  is  smaller  than  the  air 
temperature,  there  are  downward  (negative)  buoyancy  flux  and  heat 
flux.  This  tends  to  establish  a  less  neutral  lapse  rate  above  the 
surface  and  to  increase  the  temperature  variance.  In  this  case,  the 
small-scale  refractive  index  variations  are  stronger,  and  the  backs- 
catter   signal    level    will    be    larger. 


B.       THE        STRUCTURE        FUNCTION        AND        THE  TEMPERATURE 
STRUCTURE   PARAMETER 

1.       Description 

The    refractive     index    gradient     produced    by     the    interaction 

of      the      acoustic      wave      with      small      scale      velocity  fluctuations      is 
expressed    by    the    structure    function    Dn(r)    which    is    defined    by: 

Dn(r)    =    <[N(x)     -    N(x  +  r)]'>.  (3.2) 
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Here  N  is  the  index  of  refraction  at  two  points  in  space 
separated  by  a  vectorial  distance  r.  The  bracket  represents  an 
average  over  all  space.  Small-scale  turbulence  theory  predicts  that 
the  structure  function  value  is  determined  by  the  structure  param- 
eter   Cn^    and    the    scalar    separation    "r"    by    the    formula: 

Dn(r)    =    Cn'r^/'  (3.3) 

This  is  valid  only  if  the  magnitude  of  r  lies  within  the  inertial 
subrange    of    the    turbulence,     so    that    isotropy    is    "locally"    valid. 

Including  both  the  temperature  and  water  vapor  contribu- 
tions, the  structure  function  is  related  to  the  temperature  structure 
parameter  CT^,  to  the  humidity  structure  parameter  Cq^  and  to  the 
temperature -humidity    structure    parameter    Ctq    by: 

Cn^={19    10-'P/T')'(CT2  +  .113Ctq+3.2    IQ-'Cq^),  (3.4) 

where  P  is  the  pressure  in  millibars  and  T  the  absolute  temperature 
(Friehe,     1977). 

2 .       Surface    Scaling 

The  structure  parameters  can  be  related  to  the  measured 
meteorological  quantities  through  Monin-Obukhov  surface  layer  simi- 
larity   parameters    (Wyngaard,     et    al.  ,     1971    and    Wyngaard,     1973). 

CT^    =    T=''='Z-2/'f(0  (3.5) 

Cq^    =    A    Q=i=2Z-2/^f(C)  (3.6) 

where  f(0  is  the  empirical  function  found  by  Wyngaard,  et  al. 
and  updated  by  Davidson  et  al.  (1978).  The  quantity  "A"  is  a 
constant  approximately  equal  to  0.6  (Fairall,  et  al.  1980)  The 
temperature -humidity    structure    parameter    is    given    by: 
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Ctq    =    R    T-Q':=Z-'''A'/2f(0  (3.7) 

where  "R"  is  the  temperature-humidity  correlation  parameter 
(approximately  0.8  under  unstable  conditions).  The  value  of  R  is 
not  well  known  in  the  surface  layer.  Q=''  is  the  water  vapor  density 
scaling    parameter    (gm/m^)    and    is    related    to    q''=    by: 

Where    p    =    1.2    kg/m'    is    the    air    density. 

The  Monin-Obukhov  length  scale  was  defined  by  equation 
2.3  and  was  approximated,  neglecting  the  q''=  contribution,  by  equa- 
tion   2.4    ,    which    is    repeated    here: 


t    =    Z/L    =    ZkgT-    /    T^-  (3.8) 

The  function  f(C)  has  different  forms  for  unstable  (C  <  0)  and 
stable  (5  >  0)  conditions  and  is  a  constant  for  neutral  (C  =  0) 
conditions : 

f(0    =    4.9(1    -    70"^^'                 (5    <    0:     unstable)  (3.9) 

f(5)    =4.9                                            (^    =    0    :     neutral)  (3.10) 

f(0    =    4.9(1    +    2.4^^^')                (?    >    0       :     stable)  (3.11) 

The  problem  of  predicting  Cn^  is  therefore  reduced  to  finding 
values  of  q-,  T-  and  U-.  The  ratio  CTVCq^  shows  the  relative 
importance    of    T'''    and    q''': 

CT^    /    Cq'    =    T-V    A    Q-^  (3.12) 

This    gives: 

CT^    /    Cq2    =    T-V    0.6     (1.2    q-)'=    1.4    T-V    q-' 


25 


In  equation  3.4  there  is  a  factor  equal  to  3.2  10"^  between  the 
contribution  of  the  CT^  and  the  Cq^  terms.  Therefore  the  contribu- 
tion   of    Cq^    may    be    significant    (more    than    5%)    if 

3.2    lO-'CqV    CT^    >    5%. 
This    leads    to: 

1.4    T-V   ^1===^    =    CTV    Cq^    <    3.2    10-V5% 
Thus    the    final    condition    is    : 

T=V    CC--    <    0.2°K/(g/kg). 

For  a  10  %  contribution,  the  ratio  should  be  less  than 
0.1.  This  means  that  q'''  should  be  greater  than  10  times  the  value 
of  T''=  if  its  influence  is  significant.  It  is  known  already  that  we  do 
satisfy  the  necessary  condition  for  this  ratio,  and  the  value  of  Cq''= 
can  be  neglected  in  equation  3.4.  The  same  result  will  be  found 
for  the  contribution  of  Ctq,  even  with  the  imprecision  concerning 
"R". 

The  next  step  is  to  use  equations  3.9  through  3.11  in 
equation  3.5  to  'get  CT^  in  the  surface  layer  for  every  kind  of 
stability: 

CT^=    4.9T-^Z-^/^(l-7C)"^/^  K    <    0  (3.13) 

CT'=    4.9T-2Z-='^'  C    =    0  (3.14) 

CT'=    4.9T=i=2Z-2/'(l+2.44^^')  K    >    0  (3.15) 

These     three      final     equations      wUl     allow      to     determine  the 

value    of    CT^     in    the    surface    layer    with    only    the    in    situ    data.  In 

certain     circumstances      it     is     possible      to     measure      CT^      directly  by 

measuring    the    variance    of    the    temperature    at    different    points    in  the 
atmosphere,     but    this    was    not    done    for    the    data    for    this    study. 
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Sodar   Data    Reduction 


The    Data    Set 


The  sodar  data,  which  include  sodar  characteristics 
and  measurements,  were  recorded  on  an  8-inch  floppy  disk  and 
subsequently  transferred  to  9-track  magnetic  tapes  wiiich  can  be 
used      with      the      NPS      IBM      3033      mainframe      computer.  One      tape 

contains  the  data  for  the  period  going  from  11  July  at  2150Z  to  12 
July  at  llOOZ  and  another  tape  the  remaining  data  from  13  July  at 
0230Z  to  17  July  at  1130Z.  Unfortunately,  there  are  some  gaps  in 
the    data    set    due    to    technical    problems    during    the    transfer. 

(1)  Time.  The  date  and  time  are  first  obtained. 
Some  profiles  are  missing  because  of  computer  "interfacing" 
difficulties,  and  the  interval  separating  two  profiles  can  vary  a 
little  during  each  day.  This  aspect  does  not  change  the  reliability 
of  the  measurements  corresponding  to  each  time.  The  only  problem 
this  generates  is  in  the  averaging  process,  which  will  be  discussed 
in    section    B .  3 .  b . 

(2)  Gate  and  Altitude.  Each  profile  gives  a  v^alue 
proportional  to  backscattered  sound  intensity  at  different  heights  for 
the  same  time.  The  tape  gives  the  value  of  a  burst  which  is 
always  equal  to  4  ms  in  our  case.  It  also  provides  the  number  of 
bursts  per  gate,  which  allows  computation  of  the  altitude  for  each 
gate.  This  interval  was  constant  and  corresponded  to  7  bursts/gate 
for  the  first  tape  (11  and  12  July),  which  gives  an  interval  of  4.6 
meters  as  shown  earlier.  For  the  second  tape  it  varies  from  2  to 
6  bursts/gate.  This  gives  an  interval  between  1.3  and  4  meters 
from  equation  3.1.  The  maximum  heights  detectable  for  these  cases 
were,  respectively,  290  and  820  meters.  To  simplify  the  data  set 
of  the  second  tape,  it  was  necessary  to  make  a  transformation  in 
the  signal  intensity  for  each  profile  and  each  level  and  to  normalize 
the  interval  everywhere.  A  5  bursts/gate  interval  was  chosen 
because     it     was     the     value     encountered     most     of     the     time      (and     the 
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closest  to  the  other  numbers),  so  that  the  mathematical 
transformation    causes    as    little    imprecision    as    possible. 

(3)       Noise       and        Signal. .  A        background       noise 

estimate  is  obtained  during  each  profile  and  is  subtracted  from  the 
backscattered  intensity,  after  some  corrections  which  depend  on  gate 
altitude.  These  final  power  values  become  the  usable  data  for  each 
level  after  an  averaging  process  is  executed  by  the  system  for  each 
profile    and    each    gate. 

Additional  noise  can  be  caused  by  aircraft  or 
helicopter  operations.  Because  the  gain  of  the  system  is  rather 
high,  it  is  very  sensitive  to  any  external  noise  sources.  The  air 
flow  around  the  acoustic  enclosure  can  also  produce  noise.  This  was 
most    likely    when    the    wind    speed    was    higher    than    9    m/s. 

When  the  system  detected  noise  in  excess  of  an 
adjustable  threshold,  it  simply  assigned  an  arbitrary  value  to  the 
backscattered  signal  intensity  of  I  =  32767,  whereas  usual  values 
were    between    200    and    600. 

b.        The    Averaging    Processes 

Because  of  the  sometimes  important  fluctuations  of  the 
received  signal,  it  was  necessary  to  smooth  the  response  with  an 
averaging  technique.  For  this  process  as  well  as  for  the  plot  of 
the  boundary  layer,  the  programs  were  adapted  from  those  of  Lt. 
Mohn  (1985).  The  first  step  of  the  analysis  was  to  read  the 
profiles  by  groups  of  eight  and  make  an  average  for  each  gate. 
Then  we  have  a  value  of  the  backscatter  signal  ranging  from  every 
7  to  10  minutes.  This  range  is  variable  for  a  given  day  and 
increases  until  10-15  minutes  for  the  last  two  days.  As  discussed 
earlier,  the  averaging  interval  is  not  constant.  On  a  particular  day, 
there  is  not  much  discrepancy  from  an  arbitrary  constant  interval. 
However,  there  are  two  periods  which  will  be  the  object  of  some 
comments  in  the  time-height  cross-section  of  the  boundary  layer's 
section. 
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C.   REQUIREMENTS    FOR    SODAR    CALIBRATION 

The  method  to  retrieve  CT^  implies  the  calibration  of  the 
system    in    order    to    obtain    a    relation    of    the    form: 

CT^    =    f(I) 

As  CT^  and  I  should  be  proportional,  it  is  expected  to  find  a 
linear  function  of  the  first  order  for  f(I).  To  find  it,  CT^  is 
computed  in  the  surface  layer  using  equations  3.13  to  3.15  .  This 
wUl  be  done  for  the  first  six  levels,  whose  height  goes  from  25  m 
to  41.5  m  with  a  separation  of  3.3  m  given  by  equation  3.1  , 
when  there  are  5  bursts  per  gate.  This  was  not  the  case  for  the 
first  tape,  but  after  elimination  of  data  which  do  not  fit  for  this 
computation,  only  one  pair  of  data  remained  whose  value  of  CT^ 
comes    from    this    tape. 

In  the  next  step  a  program  plotting  CT^  versus  I  is  used  for 
the  six  different  heights.  To  get  good  data  for  each  set,  the  CT^ 
values  are  averaged  for  about  every  half  hour.  As  the  values  of 
the  in  situ  measurements  are  already  averaged  every  10  minutes,  3 
values  are  averaged,  keeping  the  time  of  the  middle  value.  Then 
the  values  of  the  backscattered  signal  are  read  15  minutes  before 
and    after    this    time    and    averaged    to    complete    the    set    of    two    values 

(CT^  I). 

The  quality  of  these  pairs  must  be  evaluated  before  generating 
a    CT^=    f(I)    relation,     following    some    criterion: 

•  The  height  of  the  level  considered  in  the  surface  layer  must 
be  less  than  1/10  of  the  height  of  the  inversion  layer  for  the 
equations  giving  CT^  in  the  surface  layer  to  be  still  valid.  As  will 
seen  in  figure  4.4b,  there  was  an  interesting  evolution  of  the 
inversion  layer  on  13  July.  But  this  does  not  help  very  much  to 
get  favorable  sets  of  data  because  the  inversion  layer  height  is 
always    below    250    m    from    llOOZ    to    1700Z. 

•  For  the  same  reason,  if  a  second  mixed  layer  forms  from  the 
ground     and     so     is    very     low,     or     if     we     have    a     stronger     signal     at 
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intermediate  level,  the  6  levels  considered  will  be  probably  higher 
than  1/10  this  new  height.  Therefore,  these  data  which  correspond 
to    a    level    which    is    above    the    surface    layer    must    be    rejected. 

For  example,  on  13  July  from  0300Z  to  llOOZ  it  wHl  be  seen  in 
figure  4.4b  ,  a  second  inversion  layer  whose  height  varies  from  the 
ground  to  230  meters.  The  display  terminal  shows  also  the  same 
kind  of  phenomenon  on  14  July  at  0330Z  when  a  second  inversion 
layer  begins  to  rise  up  from  the  surface  layer,  reaching  a  height 
of  about  150  meters  at  0430Z  and  200  meters  at  0530Z.  Then  it  is 
completely  mixed  with  the  surface  layer  at  0600Z.  Also  on  16  July, 
a  new  layer  forms  at  0300Z  at  a  height  of  about  250  meters.  Then 
it    lowers    to    100    meters    at    0415Z    before    joining    the    surface    layer. 

During  all  these  periods,  the  corresponding  sets  of  values  must 
be  eliminated  when  the  height  of  the  inversion  layer  is  too  low 
(below  250  meters)  and  cannot  fit  with  the  1/10  rule's  requirement 
because    the    sodar    cannot    make    any    measurement    below    25    m. 

Another  source  of  problems  is  the  fact  that  there  are  some 
gaps  in  the  data  set  of  the  in  situ  measurements  made  from  POLAR 
QUEEN  station.  This  led  to  missing  values  for  U'''.  About  half  of 
these  values  are  missing  between  11  July  at  2300Z  until  14  July  at 
2300Z.  Also  T-  and  U-  are  not  known  from  16  July  2100Z  until 
the  end  of  the  measurement  period.  Therefore,  it  is  not  possible  to 
form  pairs  of  data  for  the  calibration  during  these  intervals  of 
time. 

In  the  backscatter  signal  data  there  are  also  bad  data  due  to 
noise.  There  are  relatively  few  such  cases,  and  they  are  generally 
encountered  at  the  first  level.  For  instance,  on  14  July  from  1315Z 
untH  1700Z,  about  half  of  the  data  set  at  the  first  level  is 
unusable,  but  there  are  still  enough  values  left  for  the  averaging 
process,  which  is  done  about  every  half  hour  and  during  which  25 
profiles    are    read. 

Since  CT^  and  I  are  obtained  in  the  surface  layer  for  6 
levels,  it  is  possible  to  calibrate  the  system  after  having  deduced 
the    function    relating    CT^     and    I    at    these    low    levels.     After    that    the 
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value  of  the  backscattered  signal  needs  to  be  used  with  this  rela- 
tion to  get  the  value  of  CT^  everywhere  in  the  atmospheric 
boundary    layer. 


D.       COMPUTATION   OF    THE    INVERSION   HEIGHT 

It  is  possible  to  compute  the  inversion  layer  height  Zi  with  the 
data  set  without  generating  a  time-height  cross-section.  To  do  that 
it  is  only  assumed  that  the  inversion  height  occurs  when  the 
backscattered  signal  is  maximum  outside  the  surface  layer.  There 
are  some  difficulties  if  there  are  two  different  inversion  layers  or  if 
the  boundary  layer  is  stable  and  produces  a  substantial  return 
throughout. 

1.       The    Problem    of    Determining    the    Inversion    Layer 

The  problem  of  multiple  inversions  is  not  too  disturbing 
because  a  second  inversion  layer  is  generally  weaker  and  is  not 
accounted    for    in    the    algorithm. 

The  second  problem  is  more  troublesome.  The  strength  of 
the  signal  just  above  the  surface  may  be  stronger  than  the  inter- 
esting signal.  In  that  case,  it  is  necessary  to  read  the  data  not 
at  the  first  level,  but  at  a  higher  one,  where  the  strong  surface 
layer  return  is  not  included.  Sometimes  this  extension  from  the 
surface  reaches  higher  than  the  inversion  layer  itself.  Therefore, 
the  first  level  to  be  considered  must  be  determined,  depending  on 
the  particular  configuration  for  a  given  day  or  part  of  day.  If  not, 
the  result  will  indicate  either  the  top  of  the  surface  layer  or  the 
inversion  layer,  or  even  a  height  above  the  inversion  if  the  first 
level    considered    is    above    the    inversion. 

The  program  must  be  run  a  few  times  to  determine  the 
first  appropriate  level  to  be  read  with  regard  to  a  particular 
interval     of    time.     This    first     level    will     be    adjusted     and    changed     in 
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order    to    be    always    approximately    between    the    surface    layer    and    the 
inversion    layer. 

2.  The    Inversion    Layer    Height 

Once  this  difficulty  is  solved,  the  height  of  the  inversion 
is  computed  using  two  averaging  processes.  The  first  one,  as 
before,  consists  of  averaging  the  backscattered  signal  for  each 
level,    taking    eight    profiles    over    a    time    of    about    ten    minutes. 

The  second  process  gives  the  height  of  the  stronger  signal 
and  is  executed  by  making  a  running  vertical  average  which  spans 
five  levels.  The  stronger  value  obtained  will  give  the  level  of  the 
inversion    layer. 

3.  The    Inversion    Layer    Thickness 

Another  useful  feature  which  can  be  computed  from  the 
sodar  data  is  the  inversion  layer  thickness.  This  in  done  in  the 
same  program  that  gives  the  inversion  layer  height.  The  averaging 
processes  are  the  same,  but  instead  of  searching  for  the  strongest 
vertical  signal,  we  look  for  the  weakest  signal  below  the  inversion, 
which  should  occur  where  the  structure  function  parameter  is  the 
smallest.  This  happens  generally  at  the  base  of  the  inversion  layer 
Zl.  To  find  the  value  of  the  top  of  the  inversion  layer  Zu,  the 
level  above  Zi  is  determined  where  the  averaged  value  of  the  signal 
is  half  of  its  maximum.  This  factor  1/2  is  arbitrarily  chosen,  but 
since  the  value  of  the  signal  decreases  relatively  fast  above  the 
inversion  than  below,  a  strong  change  of  the  signal  will  occur  for 
a  small  height  difference.  Therefore,  the  imprecision  in  choosing 
this    factor    should    not    be    important. 

The    final    value    Th    of    the    inversion    thickness    is    given    by: 

Th    =    Zu    -    Zl  (3.16) 


i 
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A    listing    of    this    program    with    some    comments    is    given    in    Appendix 

"A". 
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IV.    THE   RESULTS 

In  this  chapter,  we  will  look  at  the  three  different  kinds  of 
results  obtained  with  the  methods  described  previously:  the  calibra- 
tion of  the  sodar,  the  time-height  cross-section  history  for 
MIZEX-84,  and  the  estimation  of  the  potential  temperature  and  the 
refractive    index    gradient    in    the    inversion    layer. 


A.        THE    CALIBRATION    OF    THE    SODAR 

1.       Analysis    of    the    Results 

Figure     4.1     shows     CT^     versus     I     for     the     six     backscatter 
levels    closest    to    the    surface.     Specifically    these    are: 
Figure    4.1a    :     level    1    ,        25    meters. 
Figure    4.1b     :     level    2    ,        28.3    meters. 
Figure    4.1c     :     level    3    ,        31.6    meters. 
Figure    4.  Id     :     level    4    ,        34.9    meters. 
Figure    4.1e     :     level    5    ,        38.2    meters. 
Figure    4.  If    :    level    6    ,       41.5    meters. 

a.        Level    1 

The  first  plot  seems  to  be  unreliable  because  all  the 
values  of  1  are  in  the  same  range,  whatever  the  value  of  CT^. 
The  reason  may  be  related  to  effects  of  transducer  reverberation  on 
the  gates  which  are  the  closest  to  the  sodar  system.  This  is 
suggested  by  the  fact  that  a  lot  of  values  for  I  equal  to  32767  are 
encountered  at  this  level,  an  indicator  of  signal  saturation.  An 
interesting  feature  is  the  asymptotic  approach  to  1=480  that  the 
curve  appears  to  have  as  CT^  increases.  Before  making  any 
comment,     the    other    levels    shall    be    examined. 
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Figure    4.1a  CT^    versus    Backscatter    Intensity    at    Level    1 
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Figure    4.1b         CT^    versus    Backscatter    Intensity   at    Level   2. 
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Figure    4.1c         CT^    x^ersus    Backscatter    Intensity    at    Le\>el    3. 
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Figure    4.1e         CT^    versus    Backscatter    Intensity   at    Level    5. 
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Figure    4.  If         CT^    versus    Backscatter   Intensity   at    Level    6. 
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b.  Level    2 

This  level,  shown  in  figure  4.1b,  looks  rather  like 
the  former  one.  The  reverberation  stUl  appears  to  be  significant  at 
this  range,  although  for  small  values  of  CT^,  the  intensity  is  lower 
than  before.  There  is  evidence  of  a  monotonic  functional  dependance 
of  CT^  on  I  for  400<I<475  (CT^  =  10"\  1=400)  to  (CT'  =  2  10"\ 
1=475).  After  that  point  there  is,  as  before,  a  limitation  which 
looks    like    a    saturation    of    the    system    at    1=485. 

c .  Level    3 

The    lower    points     spread    a    little    more    than    previously 
and    the    starting    point    of    the    curve    seems    to    be    at    1  =  370. 
A    few    points    seem    to    be    out    of    range: 

•  There  are  five  of  them  in  the  area  delineated  by  0<CT^<10"'*  and 
408  <  I  <  422.  One  of  these  points  is  taken  on  14  July  at  1712Z, 
when  the  inversion  layer  is  about  at  320  meters  for  the  half  hour 
average.  This  value  is  at  the  limit  of  ten  times  the  height  of  level 
3  (31.6  m) ,  and  therefore  the  bulk  formula  calculation  of  CT^  may 
be  tenuous.  It  may  also  be  due  to  non-stationarity,  in  which  case 
the  fluxes  may  no  longer  be  constant  with  height  and  the  similarity 
theory  no  longer  holds.  Advection  of  heat  and  velocity  may  play  a 
similar  role,  which  invalidates  similarity  theory  by  changing  the 
balance    of    the    terms    in    the    equations    of    motion. 

On  15  July  the  other  points  span  the  time  from  0923Z 
to  1024Z  and  occur  at  1528Z  when  the  inversion  height  is  between 
310  and  320  meters.  As  before,  these  values  are  just  barely  in  the 
limit  accepted  for  the  data  set.  This  situation  may  explain  why  CT^ 
is    smaller    than    it    should    be    given    the    value    of    the    intensity. 

•  Apart  from  this  group  of  points,  there  are  also  two  points 
(CT2=0.6  10"',  1=465)  and  (CT^  =  10"\  1=483)  which  lie  weU  below 
our  fitted  curve.  These  two  consecutive  points  occur  at  0450Z  and 
0556Z     on      16     July.         Before     0415Z      there     was     a     second      and     low 
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inversion  layer,  as  stated  earlier.  After  0415Z,  there  is  no 
apparent  reason  why  there  should  be  any  cause  to  eliminate  these 
two  values.  The  disagreement  may  be  due,  as  before,  to  non- 
stationarity,  advection  or  unusual  ambient  noise  which  lead  to 
stronger  values  of  I  compared  to  real  value  of  CT^ .  The  next 
point,  which  (CT2  =  2  10" ^  1=479)  was  obtained  at  0626Z  is  in  much 
better    agreement    with    the    fitted    curve. 

d.       Level    4 

At  this  level,  the  zero  of  the  curve  seems  to  start  at 
about  1=350.  There  are  still  two  points  which  appear  to  be  out  of 
range,  (CT^=0.6  10-^  1=459  and  CT'=0.96  10-^  1=478),  probably 
for    the    same    reasons    suggested    in    the    previous    paragraph. 

There  are  three  different  domains  evident  in  this  plot. 
The  first  represents  the  points  with  low  values  of  CT^ ,  which  are 
around  a  value  of  1=350.  Some  of  these  points  have  very  low  values 
of  the  backscattered  signal  which  occured  during  the  first  two  days 
of  the  study  period.  There  is  a  lot  of  scatter  among  the  points 
in  this  area,  probably  because  condition  were  near  neutral,  with 
small  values  of  Ta-Ts  (difference  between  air  temperature  and 
surface  temperature) .  So  the  relative  effect  of  instrumental  error  on 
Ta-Ts  may  be  very  large.  This  scatter  may  also  be  due  to  the  fact 
that,  even  for  low  values  of  CT^,  the  system  gives  too-high  values 
of    I    because    of    the    background    and    system    noise. 

The  second  area  contains  the  seven  points  around 
CT^=0.5  10"^  and  1=400.  They  all  correspond  to  the  same  period, 
from  OOOOZ  to  0320Z  on  16  July,  and  are  described  satisfactorily  by 
a  straight  line  segment.  The  third  area  includes  points  where  CT^ 
>  1.5  10"^  and  I  >  460.  Again,  all  of  these  points  occur  during 
the  period  from  0615Z  to  2040Z  on  16  July.  In  this  case,  the 
points  may  be  described  by  a  straight  line  from 

(CT'  =  1.8    10-\     1=470)       to       (CT'=4.5    10-^     1=490).  Therefore,       the 

overall      set     of      points      can     be      described      by      three     different      line 
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segments    although    the    first    one    will    not    be    very    reliable    because    of 
scattering. 

e .  Level    5 

Starting  at  this  level,  a  lot  of  profiles  whose  inver- 
sion layer  height  is  less  than  380  meters  must  be  rejected.  These 
cases  are  essentially  in  area  one  where  we  got  low  values  of  CT^ 
in  level  4.  So  little  information  is  lost  and  most  of  the  interesting 
profiles    of    level    4    are    still    there    in    that    case. 

It  can  be  noted  that  the  slope  of  the  fit  looks  less 
curved    than    for    previous    levels. 

f.  Level    6. 

This     level     gives  results     similar     to     level     5,     although 

we    do    not    have    a    well-defined  straight    line    in    the    second    area.     But 

the     general     tendency    remains  the     same     and    shows     good     agreement 
with    the    previous    plots . 

2.       Calibration    of    the    System 

To  obtain  a  reliable  calibration  curve,  the  plots  that  are 
likely  to  be  free  of  reverberation  effects  will  be  kept.  These  corre- 
spond to  the  last  three  levels  (4,  5  and  6).  Here  the  points 
corresponding  to  the  profiles  taken  at  0450Z  and  0556Z  on  16  July 
are  eliminated  for  the  reasons  described  earlier.  The  combination  of 
these    three    plots    is    shown    in    figure    4.2. 

The  three  sets  of  data  superimpose  quite  well  but  do  not 
fit  on  a  straight  line  as  expected  if  there  is  proportionality  between 
CT^  and  1.  However,  it  is  useful  to  plot  the  same  data  in  loga- 
rithmic coordinates,  log(CT^)  versus  log(I),  to  investigate  whether  a 
power  relation  provides  a  better  fit  for  the  relation  between  CT^ 
and     I.         The      resulting     plot     is     shown     in     figure      4.3.         The     two 
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Figure    4.2         CT^    \^ersus    Backscatter    Intensity 
at    Levels    4,    5    and    6. 
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Figure    4.3         Log(CT^)    versus    Log(I)    at    Levels    4,    5    and    6. 
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regions      with     higher      backscattered      signal     fit  rather      well     on      a 

straight     line.      The     least-squares     approximation  gives     the  following 
equation : 

LogCCT^)    =    8.63    Log(I)    -    25.66  (4.1) 

For    lower    values    of    I,    there    is    much    scatter    in    the    data.  This    is 

more    apparent    on    this    figure.        Therefore,     the  calibration  given    by 

equation    4.1    should    be    used    carefully    for    I    less  than    380. 


3.       Error   Analysis 

To    find    the    relative    error    on    CT^     the    value    of    the    vari- 
ance   o^    is    first    computed    using    the    approximation: 

o^    =    1/N    Z(log(CT2    -    8.631og(l)     +    25. 66)^ 
where    N    is    the    total    number    of    points     considered.        o^     is    found    to 
be     equal     to     0.01.         The     relative     standard     error     on     log(CT^)      is 
given    by: 

d(log(CT^))/log(CT^)    =    aVlog(CTM 
whose    left    hand    side    is    also    equal    to: 

d(ln(CT')/ln(CT')    =    d(CT')/(CT'    InCCT^)) 
These    expressions    give: 

d(CT')/CT2    =    0^    ln(CT')/log(CT') 
Finally,     the    relative    uncertainty    on    CT^    is: 

d(CT2)/CT^    =    o^    log(lO)    =    23% 


B.       TIME-HEIGHT    CROSS-SECTIONS    OF    THE   INVERSION   LAYER 

There  are  several  gaps  in  the  time  history  of  the  sodar  data, 
so  the  figures  showing  the  boundary  layer  and  its  inversion 
obtained    from    sodar    data    are    decomposed    as    follows: 

Figure    4.4a    :     11    July    2150Z    to    12    July    llOOZ. 

Figure    4.4b     :     13    July    0230Z    to    2400Z. 
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Figure  4.4c  :  14  July    OOOOZ  to  2400Z. 

Figure  4.4d  :  15  July    OOOOZ  to  2400Z. 

Figure  4.4e  :  16  July    OOOOZ  to  2400Z. 

Figure  4.4f  :  17  July    OOOOZ  to  1130Z. 

The     data    are    placed    in    a     two    dimensional    array,     one    carrying 

the  time,  the  other  the  altitude  (or  the  gate).  The  lowest  value  of 
the  backscatter  intensity  to  be  plotted  is  determined  by  a  value 
"BSMIN".  This  is  chosen  for  each  plot  such  that  the  isopleths  are 
so  crowded,  and  the  surface  layer  and  the  interfacial  layer  are 
most  apparent.  The  lower  the  strength  of  the  inversion  layer,  the 
lower  the  value  of  BSMIN.  An  incremental  value  is  also  chosen  so 
that  the  contour  lines  are  not  too  tightly  spaced.  The  plot  gives 
essentially  a  contour  of  the  backscattered  signal,  but  the  values  of 
CT^    according    to    the    previous    calibration    have    been    added. 

The  resulting  plots  are  first  examined  and  their  behavior  are 
discussed  as  it  relates  to  environmental  conditions.  Before  doing 
this,  it  must  be  recalled  that  the  summer  "night"  in  these  latitude 
is  much  like  day.  In  fact,  the  downward  radiation  during  the 
night    is    roughly    equal    to    2/3    of    the    value    during    daytime. 

1.       11    and    12    July 

This  period,  which  includes  the  night  between  days  11  and 
12,  shows  a  smooth  and  continuous  decrease  in  the  inversion  layer. 
This  is  associated  with  a  more  stable  layer  which  characterized  the 
night  time  and  a  relatively  high  pressure  (1013  mb)  on  12  July 
between  1200Z  and  2400Z  which  caused  subsidence  and  tended  to 
lower  the  inversion  altitude.  The  low  value  of  the  chosen  BSMIX 
(290)  indicates  that  the  strength  of  the  inversion  is  not  very 
great.  The  background  noise  for  the  period  is  about  230,  not  very 
far    below    290. 


47 


M 

LO 

O 

rv 

O 

— * 

1 

t— 1 

!— 

_1 

2 

=3 

LJ 

— ! 

IT 

UJ 

CM 

ce 

CJ 

z 

O 

f-i 

M 

O 

O 

OJ 

(M 

_] 

=3 

— ) 

,^^ 

r—t 

o 

CD 

• 

OO 

■^ 

CD 

i 

X 

3 

LJ 

<— 4 

M 

2= 

on 

^ 

CD 

» 

^ 

o 

Q 

— ^ 

-^ 

LO 

in 

o 

ro 

CJ 

^ 

^ 

CD 
O 

in 

'-' 

•  0  CZ 


300- 
0 


bo 


0^ 


^^s?. 


!'> 


N-V 


0  oo^^^y:^:xr>  •  ^ 


0      s^S^J 
0  0  ♦♦        c> 


^^0 


"r^'l 


CD 


o{)\bocu3:>- 


♦     ♦     ♦  0 


r-  a,' 


^  o 


300 -^^ V/^V/^''~^ 

-^  -inn  / 


♦  0 


SS9 


SSS 


S2Z 


Figure    4.4a 


Atmospheric    Boundary    Layer    Structure 
on    11    and    12    of   July. 
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Figure    4 .  4b 


Atmospheric    Boundary    Layer    Structure 
on    13    of   July. 
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Figure    4 .  4c 


Atmospheric    Boundary    Layer    Structure 
on    14    of   July. 
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Figure    4 . 4d 


Atmospheric    Boundary    Layer    Structure 
on    15    of    July. 
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Atmospheric    Boundary    Layer    Structure 
on    16    of   July. 
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Figure   4 .  4f 


Atmospheric    Boundary    Layer    Structure 
on    17    of   July. 
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2.  13    July 

In  the  second  half  of  the  night  we  see,  as  before,  a 
decrease  in  the  inversion  layer.  The  stability  which  tends  to  occur 
contributes  to  the  development  of  a  period  of  intermittent  strong 
backscatter  at  intermediate  levels  which  rose  and  mixed  with  the 
higher    one    after    sunrise. 

Until  noon,  the  wind  speed  was  low  and  U=''  was  about  10 
cm/s.  Then  U'''  increased  to  40  cm/s  and  remained  between  30  and 
40  cm/s  until  1800Z.  This  indicates  a  stronger  mechanical  mixing, 
which  raises  the  inversion  layer  through  entrainment.  Of  course  this 
increase  may  also  occur  as  a  response  to  the  incoming  solar  radia- 
tion, which  generated  more  unstable  conditions  and  therefore 
buoyant    mixing    during    the    day. 

Dorman  (1985)  suggests  that  gravity  waves  may  contribute 
to  the  kind  of  movements  of  the  inversion  that  we  can  see  during 
this  day.  This,  as  well  as  advection,  can  cause  the  deepening  and 
then  the  raising  of  the  inversion  layer.  It  is  difficult  here  to  have 
a  quantitative  appreciation  of  these  effects  because  of  the  lack  of 
supporting  environmental  data.  A  higher  value  of  BSMIN  (320) 
indicates    the    strength    of    the    inversion    layer    for    this    day. 

3.  14    July 

There  were  no  data  between  0503Z  and  0911Z.  This  caused 
a  jump  in  the  time  scale  for  the  plot  of  the  inversion  layer.  This 
plot  has  been  smoothed  by  the  averaging  process.  In  fact  the  true 
inversion  height  goes  from  500  meters  at  0500Z  to  390  meters  at 
0900Z. 

There  was  another  gap  between  2028Z  and  2203Z,  but  the 
inversion  stayed  nearly  at  the  same  height  during  this  interval. 
This  is  why  the  horizontal  scale  of  figure  4.4c  does  not  show  an 
interruption    for    this    short    period. 

As  before,  a  stronger  backscatter  region  developed  at 
intermediate    level     from    0330Z    until    0600Z.     This    is    not    easy     to    see 
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on  the  plot  because  of  the  lack  of  data,  but  in  this  case  it  rises 
up  to  180  meters  then  disappears  with  the  early  morning.  This  may 
be  caused  by  large-scale  processes  such  as  waves  or  advection. 
This    day    is    also    characterized    by    a    weaker    inversion    layer. 

4.  15    July 

There  are  no  data  for  a  one  hour  period  between  1620Z 
and  1727Z.  The  inversion  layer  went  from  265  meters  to  245  meters, 
which  is  a  decrease  of  less  than  10%.  Therefore  there  is  no  need 
for  an  interruption  in  the  horizontal  scale  of  figure  4.4e,  although 
the    reader    should    be    aware    of    this    gap    period. 

There  was  not  much  change  during  this  24  hour  period. 
The  signal  in  the  inversion  layer  has  a  value  greater  than  450. 
This  and  the  13  July  data  were  the  two  periods  where  the  inver- 
sion   layer    will    have    probably    the    strongest    temperature    gradients. 

5.  16    July 

The  heat  flux  was  negative  (downward)  during  this  entire 
day,  and  conditions  were  therefore  stable.  The  temperature  on  the 
ice  mast  at  5.5  meters  went  from  0.2°C  at  OOOOZ  to  1°C  at  noon 
and  up  to  2°C  at  2400Z.  The  temperature  just  above  the  ground 
was  around  1°C  from  lOOOZ  to  1700Z  and  continued  to  increase  to 
1.7°C  at  2400Z.  These  values  of  the  ambient  temperature  were 
slightly    higher    than    the    other    days. 

This  increase  of  the  net  heat  in  the  boundary  layer 
should  be  responsible  for  the  continuous  growth  of  the  inversion 
layer  during  all  the  day.  It  is  not  easy  to  see  what  role  the 
subsidence  played  because  we  were  in  a  period  of  high  pressure 
(1016  mb),  with  subsidence  expected  to  lower  the  inversion  layer 
height.  Obviously,  this  was  not  the  case  here.  The  wind  speed 
remained  between  4  and  6  m/s  until  noon  and  then  increased  and 
stayed     between     6     and     9     m/s.      This     wind     caused     more     mechanical 
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mixing     and     may     also     be     responsible     for     the     increasing     inversion 
layer    height. 

6.       17    July 

The  wind  continued  to  be  relatively  strong  during  this 
day.  It  stayed  between  7  and  10  m/s  and  is  probably  the  reason 
why  the  inversion  layer  continues  to  behave  as  it  did  during  the 
day  before.  There  are  no  in  situ  data  of  U'''  and  T'"  for  this 
period  to  get  the  value  of  the  heat  flux.  The  temperature  just 
above  the  ice  and  at  5.5  meters  are  very  close,  so  conditions  must 
have  been  near-neutral.  This  can  explain  the  weak  inversion  layer 
which  tends  to  vanish  and  the  greater  interval  in  the  surface  layer 
between  the  contour  lines  than  on  the  previous  days.  This  was  a 
time  of  formation  of  an  area  of  stronger  backscattered  signals  at 
intermediate  levels  around  03-0400Z  and  0900Z,  which  decayed 
subsequently. 


C.       ESTIMATION       OF        POTENTIAL       TEMPERATURE        JUMP       AND 
REFRACTIVE-INDEX    GRADIENT    AT    THE    INVERSION    LAYER 

1.       The    Potential    Temperature    Jump 

Wyngaard  and  Lemone  (1980)  derived  an  equation  which 
shows  that  it  is  possible  to  relate  the  mean  value  of  the  virtual 
temperature  structure  parameter  with  the  jump  of  potential  virtual 
temperature  at  the  inversion  layer,  the  height  of  the  inversion  and 
the    virtual    temperature    scaling    parameter: 

<CTv^>    =    0.5    Tv-    ATv    /    Zi^/'  (4.2) 

where  the  brackets  mean  average  over  the  inversion  layer  and  the 
subscript  v  indicates  virtual  temperature.  This  is  only  valid  for 
convective       cases.       However,       since      there      were       no       strong      free 
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convection  conditions  in  the  Arctic,  we  will  consider  unstable  cases 
for    values    of    Qo    greater    than    10"-'    °K    m/s. 

Wyngaard  and  Lemone  also  show  how  to  estimate  the 
potential  temperature  jump  at  the  inversion,  assuming  that  the 
surface  moisture  flux  and  the  mean  humidity  jump  can  be  neglected 
at  the  top.  These  assumptions  allow  use  of  normal  temperature 
instead    of    virtual    temperature. 

Tv'!=    =    T=!' 
Tv       =    T 

The  first  assumption  is  reasonable  because  q'''  is  very  small  in  this 
experiment,  so  the  vapor  flux  U''=q'''  will  be  negligible.  Tlie  second 
assumption  is  not  so  reasonable  because  the  value  of  the  humidity 
jump  at  the  inversion  is  not  known.  The  presence  of  fog  in  the 
area  is  an  indicator  of  moisture,  but  it  can  be  expected  that  the 
cold  temperatures  (0°C)  will  produce  a  generally  drier  atmosphere 
(due  to  lower  saturation  vapor  pressure) .  Radiosonde  data  wUl  be 
helpful  in  that  case  but  are  not  currently  available.  Therefore,  the 
deduction  of  the  inversion  layer  potential  temperature  jump  should 
be  taken  as  a  first  guess.  With  these  qualifications,  equation  4.2 
can    be    rewritten: 

AT    =    2    <CT^>    Zi^/V    T'-  (4.3) 

The  only  problem  now  is  to  compute  the  mean  of  CT^  in  tlie  in\'er- 
sion  layer.  Wyngaard  and  Lemone  recommend  assuming  a  peak-to- 
mean  ratio  for  CT^  of  2.  The  peak  value  for  CT^  is  easUy 
determined  because  we  know  Zi  and  the  corresponding  value  of  the 
backscatter  signal.  Then  equation  4.1  is  used  to  get  the  value  of 
CT^  at  the  peak  (CTp^),  which  allows  to  infer  the  potential 
temperature    jump    in    the    inversion.        Equation    4.3    becomes: 

AT    =    CTp^    Zi'^V    T-=  (4.4) 
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2.       Refractive    Index    Gradient   at    the    Inversion    Layer 

a.       Derivation    of    the    Refractive    Index    Gradient 

The  value  of  the  refractive  index  gradient  at  the 
inversion  layer  can  be  obtained  from  the  equation  for  refractivity 
(Meeks,     1982): 

N    =    77.6    P/T    +    3.73    lO'e/T^  (4.5) 

Where  P  and  e  are  the  atmospheric  pressure  and  the  water  vapor 
pressure  respectively  in  millibars  and  T  is  the  temperature  in  °K. 
Assuming  that  e  is  negligible  which  is  compatible  with  the  previous 
assumption  to  derive  AT,  and  taking  the  derivative  with  respect  to 
Z,    equation    4.5    becomes: 

dN    /    dZ    =    m    +    n    dT/dZ  (4.6) 

Where : 

m    =     (77.6/T)    dP/dZ 

n    =    -77.6/T^ 
P     at     the     surface     was     approximately     equal     to     1005     mb.      The     ratio 
dP/dZ    is    given    by    the    hydrostatic    approximation: 

dP/dZ    =    -    Pg/RT  (4.7) 

where  R  is  the  gas  constant.  For  normal  atmospheric  conditions  R 
=  287  J°K-'kg-^  Assuming  that  T  =  273°K,  we  get  dP/dZ  =  -125 
mb/km.  The  inversion  height  was  typically  400  meters,  which  gives 
P  =  960  mb  at  the  inversion  and  the  new  value  of  dP/dZ  in  the 
inversion    layer    reduces    to    -120    mb/km. 

Then    the    values    of    the    two    parameters    m    and    n    are: 

m    =    -34    km"^ 
n    =    -0.001    °K-^ 
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Equations  3.16,  4.4  and  4.6  give  the  value  of  dN/dZ 
at    the    inversion: 

dN/dZ    =    -34    -    (CTp^    Zi^  '  V    T-    Th)  (4.8) 

Where    Zi    and    Th    are    in    meters    and    dN/dZ    in    km~\ 

b.  Refractive    Index    Gradient    and    Refraction 

Hitney       et       al.  (1981)        developed       the       Integrated 

Refractive  Effects  Prediction  System  (IREPS),  which  is  used  to 
determine  the  radar  coverage  in  an  area  as  a  function  of  the 
refractive  condition  determined  by  atmospheric  parameters  as  pres- 
sure temperature  and  humidity.  IREPS  classification  defines  the 
following    conditions    of    refraction    with    respect    to    dN/dZ: 

Subrefraction:  0    <    dN/dZ 

Normal:  -79    <    dN/dZ    <    0 

Super    refraction:  -157    <    dN/dZ    <    -79 

Trapping:  dN/dZ    <    -157 

The  modified  index  of  refraction  M  is  related  to  N  by 
the    formula : 

M    =    N    +    157    Z  (4.9) 

Where  Z  is  the  height  in  km.  Therefore  we  will  have  trapping 
when    dM/dZ    <0    or    AM    <0    for    a    positive    dZ. 

c.  Computed    Refractive    Conditions 

The  values  of  AT,  dN/dZ,  AM  and  some  variables  used 
to  retrieved  these  values  are  given  in  appendix  "B".  figure  4.5 
show  the  profile  of  AM  versus  time  in  the  inversion  layer.  The 
trapping  conditions  happen  for  the  points  which  are  below  the  hori- 
zontal   line     (AM=0)  .        Some    particular    points    have    been    emphasized    by 
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a  marker.  They  depicted  the  largest  10  %  of  Qo  values.  Under 
these  conditions,  the  computation  should  be  more  accurate  because 
these    values    are    closest    to    the    free    convection    case. 

The  refraction  is  generally  normal  until  1800Z  on  14 
July.  After  this  period  and  until  2200Z  on  15  July  there  are  much 
stronger  values  for  dN/dZ  as  well  as  for  AT.  But  there  are  some 
important  fluctuations  in  the  range  of  this  data  and  the  values 
obtained  for  AT  look  sometimes  very  high  for  the  corresponding 
thickness  of  the  layer  which  stays  between  80  and  150  meters. 
These  fluctuations  may  be  due  to  the  different  errors  and  appro.xi- 
mations  introduced  in  the  computation.  The  more  sensitive  errors 
are  recalled:  error  on  Zl  and  Zu  which  determine  the  thickness, 
error  on  calibration  of  the  sodar  system  which  is  use  to  get  CTp^ 
and  approximation  concerning  the  mean  humidity  jump  at  the  inver- 
sion and  the  1/2  factor  to  retrieve  the  mean  value  of  CT^  from  its 
peak  value  in  the  inversion  layer.  The  major  error  is  probably 
the  free  convection  assumption.  The  highest  values  of  AT  occurred 
very  often  when  Qo  is  smallest,  which  means  that  the  situation  is 
less  unstable  and  that  equation  4.4  is  probably  not  valid,  which  is 
not  surprising.  All  these  features  indicate  caution  is  required  in 
interpreting  these  results.  However,  it  is  possible  to  deduce  some 
general  tendencies  during  the  covered  period.  There  should  be  more 
refractive  conditions  during  15  July  with  some  peaks  around  0200Z 
and  between  lOOOZ  and  llOOZ.  This  stronger  refraction  was  expected 
in  section  "B.4"  when  looking  at  the  time-height  cross-section, 
because  of  the  stronger  values  of  the  backscattered  signal  in  the 
inversion  layer.  On  the  other  hand  the  same  e.xpected  result  is 
not  obtained  for  the  13  July,  which  shows  only  two  periods  (1642Z 
and    1702Z)    with    dN/dZ    smaUer    than    -79    km~^ 

The  restriction  concerning  the  validity  of  the  computa- 
tion for  unstable  cases  only,  does  not  allow  to  have  any  data  after 
15    July    at    2200Z,     because    then    Qo    becomes    too    small. 
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3.       Comparison    of    some    Results    with    Radiosonde    Data 

At  the  time  this  thesis  was  written,  there  were  only  two 
radiosonde  profiles  available.  They  are  showed  in  figure  4.6  and 
were  taken  from  13  July  at  llOOz  and  14  July  at  1103Z.  The  left 
profiles  represent  temperature  and  dew  point  with  respect  to  height 
whereas    the    right    profiles    depict    the    wind    speed    and    direction, 

a.  Comparison    with    Radiosonde    of    13    July    at    llOOZ 

The  upper  and  lower  levels  of  the  inversion  measured 
on  figure  4.6  give  respectively  370  and  240  meters,  which  result  in 
a  thickness  of  130  meters.  The  corresponding  AT  is  equal  to 
3.5°K.  The  program  which  computes  Zu  and  Zl  gave  us  366  and 
273  meters  at  1057Z  and  369  and  204  meters  at  1107Z,  respectively. 
Therefore,  the  two  values  of  the  upper  level  are  in  very  good 
agreement.  For  the  lower  level  and  after  an  extrapolation  between 
1057Z  and  1107Z  to  find  the  value  at  llOOZ  (about  1/3  of  the 
interval) ,  a  value  of  230  meters  is  found  which  is  also  in  good 
agreement    with    the    radiosonde    data. 

The  values  of  AT  have  been  obtain  with  equation  4.4 
by  reading  directly  the  values  of  T'''  from  the  in  situ  data  set  and 
the  values  of  I  and  Zi  for  the  other  data  set.  Each  of  these  are 
averaged  about  every  10  minutes.  It  is  assumed  that  there  is  coin- 
cidence when  the  time  between  the  two  data  sets  is  less  than  5 
minutes.  The  value  of  AT  from  1047Z  to  1117Z  are:  3.4,  2.8,  5.6, 
2.6,  the  closest  being  equal  to  2.8  at  1057Z.  This  is  fairly  good 
agreement    with    the    3.5°K    given    by    the    radiosonde. 

b.  Comparison    with    Radiosonde    of    14    July    at    1103Z 

Proceeding  as  before,  we  look  at  the  set  of  values  surrounding 
1103Z.  The  first  one  is  at  llOOZ  (Zu=  382  m,  Zl=  336  m)  and  the 
other  one  at  1107Z  (Zu=  382  m,  Zl=  253  m) .  This  gives  at  1103Z: 
382    m    for    the    upper    level    and    300    m    for    the    lower.     The    radiosonde 
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measurement  gives  respectively  510  and  320  meters.  The  base  of  the 
inversion  seems  to  be  correct,  but  at  the  top  there  is  a  big 
difference.  A  reason  which  can  explain  this  discrepancy  may  be 
found  on  the  display  terminal  whose  shaded  plot  is  reproduced  in 
figure  4.7.  after  llOOZ,  just  above  the  marked  dark  inversion 
layer,  another  shaded  area  can  be  seen  which  may  be  a  weak 
prolongation  of  the  inversion  layer.  This  is  probably  too  weak  to 
be  detected  by  the  program,  although  the  radiosonde  responded  to 
it.  Also,  the  decay  of  turbulence  above  Zi  may  be  an  explanation 
of  this  feature.  This  less  turbulent  part  is  associated  with  a 
weaker  gradient  of  potential  temperature  than  is  the  more  turbulent 
part  around  Zi.  Of  course,  this  will  disturb  the  computation  of  AT 
and  dN/dZ.  In  fact,  the  radiosonde  gives  a  value  of  9.1°K  while 
the  program  gives  3.3°K  at  1053Z,  1.4°K  at  llOOZ  and  1107Z  which 
are    out    of    range    of    the    measured    value. 

Two  comparisons  are  not  enough  to  verify  the  method, 
but  as  the  occurrence  of  the  extended  weak  inversion  layer  is  not 
common  in  the  data  set,  we  can  expect  to  have  a  good  agreement 
between  the  computed  variables  and  the  radiosonde  data  in  most 
cases    if    we    are    in    free    convection    conditions. 
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V.  CONCLUSIONS    AND    RECOMMENDATIONS . 


The  analysis  of  one  particular  week's  data  during  the  last  days 
of  MIZEX  84  was  made  for  its  well-formed  inversion  layer.  This 
analysis  allowed  finding  a  power  law  to  calibrate  the  sodar  system 
which  was  used  to  study  the  atmospheric  boundary  layer.  The 
method  gives  only  two  periods  where  the  pairs  of  values  CT^  and  1 
were  reliable  for  this  calibration.  Both  are  during  16  July,  from 
OOOOZ  to  0300Z  and  from  0615Z  to  2040Z.  In  the  other  cases,  either 
there  were  no  in  situ  data,  or  the  inversion  layer  was  too  low,  or 
we  had  formation  of  a  second  low  inversion  layer  which  does  not 
allow  use  of  the  usual  formulas  to  compute  CT^  in  the  surface 
layer.  It  wlLL  be  interesting  in  subsequent  analyses  to  extend  the 
method  to  other  days  of  the  experiment  in  order  to  see  if  the  cali- 
bration   remains    in    agreement    with    the    relation    found    in    this    study. 

This  calibration  does  not  seem  to  hold  for  values  of  I  lower 
than  380  (or  CT^  =  4  10"'').  However,  this  is  a  minor  problem 
because  the  values  of  CT^  in  the  inversion  layer  where  the  backs- 
cattered  signal  has  a  value  generally  higher  than  400,  were  mainly 
used    in    this    study. 

The  time-height  cross-sections  give  reliable  pictures  of  the 
boundary  layer  and  especially  of  the  inversion  layer.  The  numerical 
indications  provide  some  complementary  quantitative  information  on 
the  value  of  the  temperature  structure  parameter  in  the  inversion 
layer. 

The  calibration  relation  was  then  used  to  deduce  values  of  CT^ 
to  estimate  the  potential  temperature  jump  and  the  refractive  index 
gradient  in  the  inversion  layer.  The  error  is  of  order  23  %.  Free 
convection  formulations  indicate  that  there  should  be  trapping  condi- 
tions on  15  July  around  0200Z  and  10-llOOZ.  However,  since  the 
computations    are    valid    only    for     free    convection,     these    results    must 
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be  considered  with  caution  because  Qo  was  never  very  large,  about 
10"^  °K/(m/s)  for  the  trapping  condition  period.  The  highest  value 
of  Qo  was  around  2  10"  ^  °K/(m/s)  on  13  July  between  1400Z  and 
1500Z,  but  the  refractive  index  gradient  indicated  a  normal  propa- 
gation   condition    in    the    inversion. 

The  strength  of  the  backscattered  signal  in  the  inversion 
should  have  lead  to  a  trapping  condition,  especially  on  the  13th 
and  15th  of  July,  according  to  the  time-height  cross-sections.  This 
was  confirmed  for  the  15  July,  with  the  computation  of  the  refrac- 
tive index  gradient.  It  was  apparently  not  true  on  13  July, 
although  Qo  was  higher  on  that  day  but  probably  not  high  enough 
to    be    a    free    convection    case    as    required    by    the    theory. 

The  comparison  of  free  convection-determined  AT's  with  some 
radiosonde  data  was  quite  good  in  one  case  and  erroneous  in  the 
other  because  of  an  higher  inversion  layer  top  than  expected  from 
sodar  data.  To  make  a  better  comparative  analysis,  it  will  be 
necessary  to  obtain  more  radiosonde  data.  This  should  be  a 
productive    avenue    for    future    analysis. 

Another  source  of  error  of  which  we  were  aware  too  late  to 
consider  was  that  the  value  of  I  given  by  the  sodar  data  is  in 
fact  the  summation  of  the  backscattered  power  and  a  spreading  loss 
correction    term    AI    which    is    equal    to: 

AI    =    ln(t/4)^    -    50 

where  t  is  the  time  in  ms  for  a  pulse  to  propagate  from  the 
antenna  to  the  middle  of  a  gate.  Therefore  AI  is  a  function  of  alti- 
tude and  goes  from  -44  (for  z  =  25m)  to  about  -38  when  z  = 
500m.  This  represents  approximately  a  10  %  error  in  the  surface 
and  in  the  inversion  layer.  This  should  not  cause  a  significant 
change  in  the  overall  results,  although  it  should  be  incorporated 
into    further    studies    using    this    sodar    system. 

This  kind  of  sodar  study  could  be  extended  to  other  areas 
with  stronger  inversions  and  surface  heat  fluxes  and  could  be  a 
good  means,  even  for  an  isolated  platform,  to  detect  ducting  condi- 
tions   when    a    radiosonde    is    not    available. 
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APPENDIX    A 
COMPUTATION    OF    THE   INVERSION    LAYER    HEIGHT    USING    SODAR 

DATA. 


//ZIBB    JOB    (0985,0262)    'ROUGE',  CLASSIC 

//-MAIN    ORG=NPGVM1.0teP 

//    EXEC    FORTXCG 

//FORT.SYSIN    DD    - 

C 

This    program    reads    sodar    data    from    MASS    STORAGE 
Date,     time  backscatter    signal.     Then    it    averages 

the    signal    (INT)    among    a    number    "NGATE"    of    values 
keeping    the    time    of    the    middle    value. 


C 
C 


Then    it    looks    in    one    profile    for    the    stronger  * 
value    of    backscatter    signal    (average    over      NZ" 
vertical    values)    which    Determines    the    value    of    the 
inversion    layer    height    "ZI"     (if    signal    is 
greater    than    a    certain    value    "BMIN").     It    also 

jht    of 


deter 


mines    the    lower    and    upper    heigl 


the    layer. 


VARIABLES: 

NGATE    =    No.     of    profiles    to    average. 

NVERT    =    No.     of    gates    to    consider. 

IFL  =    First    gate    ti    read    backscatter    signal    to 

compute    "ZI".     This    allow    to    eliminate    strong 
values    encountered    in    the    surface    layer. 

BBB  =    Minimum    backscatter    signal    level    to    consider. 

NZ  =    Number    of    vertical    value    to    average. 

BRUIT    =    Average    value    of    the    signal    over    NZ    levelS . 

BMIN       =    Maximum    average    value    of    the    signal    over    NZ. 

BIBI,BMAX    :     BIBI    Determines    the    maximum    backscatter 
signal    level    (BMAX)    in    the    inversion    layer. 
=    Minimum    average    value    of    the    signal    over    NZ. 


BMOU 

BMAR 

ZLl 

DD 

OUT 

V 

ZI 

ZL 

SIG 

ZU 

SIGU 

MAI 

THI 


=    Corresponding    height    of    BMAR. 

=    Number    of    bursts    per    gate. 

=    Out    of    range    variable. 

=    Speed    of    sound    in    the    boundary    layer 
=    Inversion    layer    height. 

=    Height    of    the    lower    boundary    of    inv.     layer. 
=    Corresponding    signal    of    ZL. 

=    Height   of    the    upper    boundary    of    inv.     layer, 
=    Corresponding    signal    of    ZU. 
=    Level    where    we    oegin    to    look    for    ZU. 
=    ZU    -    ZL:     Thickness    of    the    inversion    layer 


INT    (NGATE,  256)    =    Value    of    the    signal    read    in    MSS. 
DIMENSION    INT    (8,256),     BAC(280),     ZL1(280),     BMAR(280) 


NZ 

=    5 

NVERT 

=    160 

IFL 

=    85 

NGATE 

=    8 

BBB 

=    150. 

OUT 

=    -999. 

V 

=    331.0 

92 


-■"•=       Write    the    label 
WRITE    (6,     92) 
FORMAT(lX,58('-')/lX,'|DATE 


TIME  I        SIGNAL 


ZI    I', 
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c 
c 
c 
c 
c 


•signal  zlow  I  signal  ZUP| 


|71X,58('-')) 


111 


c 
c 
c 
c 
c 
c 
c 


112 

10 

13 


107 


C- 

C 

C 

C 

C 

C- 


113 


There    are    two    tapes    to    read    in    MSS 
UNIT    16:     11  12    July,     FLAG    =    1 

UNIT    11:     13    TO    17    July,     FLAG    =    2 

FLAG    =    1 

DO    42    L=    1,     2222 

DO    13    1=    1_,    NGATE 

IF    (FILaG.EQ.2)    go    to    112 
READ    (16^10,END=998)(INT(I,J),J=1, 
GO   TO    i:^ 
READ    fll,10,END=999)(INT(I,J),J=l, 
FORMAT    (       2(12815)        ) 
CONTINUE 


256    ) 
256    ) 


We    need    to    make    some    adjustements    to    the    first 
level    to    consider    to    avoid    the    confusion    with 
the    surface    layer.     The    number    below    are    choosen 
after    some    run    and    looking    at    the    display 
terminal    for    each    period. 


IF  (FLAG.EQ.l) 
IF(INT(4,4 
IF(INT(4,4 
IF(1NT(4,4 
1F(INT(4,4 
IFaNT(4,4 


GO    TO    113 

.  EQ .  13 .  AND .  INT(4,  5) .  GT.  10 
.GE.13.AND.INT(4,5).GT.21 
.  EQ .  14 .  AND .  INT  (4 ,  5 ) .  GE .  10 
.GE.15)  IFL  =  60 
.GE.16)  IFL  =  110 
IF(INT(4,4).GE.16.AND.INT(4,5).EQ.5) 
IF  INT(4,4).EQ.17.AND.INT(4,5).GT.5) 


IFL=  40 
IFL=  70 
IFL=  80 


1FL=  135 
IFL=  120 


Make  the  average   (BAG)  of  the  signal  (INT) 


DO    38    J    =    1,NVERT 
TOT    =    0.0 
N    =    0 

DO    22    I=l_,    NGATE 

K=    J^    +    IFL    -    1 


■''■''       Disregard    too    big    values    in    the    averaging    process 
=:-:=       "X"    means    we    compute    average    for    the    "l\tn"    leve 


C 
C 
C  =•'-'•=       'K"    means    we    comp 

c  r. 

c 

IF    (INT(I,K).GT.700)    GO    TO    22 

N    =    N    +    1 

TOT=    TOT    +    FLOAT  (    INT(I,K)     ) 
22  CONTINUE 

C 

IF    (N.NE.O)    GO    TO    31 
BAC(K)    =    OUT 
GO    TO    38 
31  BACfK)    =    TOT    /    FLOAT(N) 

38  CONTINUE 

Q 

C 

C  =!-       Compute    "    ZI    " 

C  ==^=========== 

c 

Q 

DO    41    JK    =    1,     280 

ZLl(JK)    =    0.0 

BMAR(JK)    =    0.0 
41  CONTINUE 

BMIN  =    BBB 

BMOU  =    500. 

ZL  =0. 
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Look    for    the    lowest    value    of    the    signal 
corresponding    to    a    height    "ZLl"    to    determine    the 
lowest    boundary    "ZL"    of    the    inversion    layer. 


NNVERT    =    NVERT    -    4 

IF    (FLAG.EQ.l)    DD    =    7. 
IF    (FLAG.  £0.2)    DD    =    5. 
DO    55    1    =    1^    NNVERT 
BIBI    =    BfeB 
AVG    =    0. 
N    =.   0 
C 
C  =;-:=       Average    signal    over    "NZ"    levels. 

C 

DO    50    J    =    1,    NZ 

K=    IFi    +    J    +    I    -    2 
IF    (BAC(K).EQ.OUT    )    GO    TO    50 
IF    (BAC(K).GE.BIBI)    BIBI    =    BAC(K) 
N    =    N    +    1 
AVG=    AVG    +    BAC(K) 
50  CONTINUE 

C 

IF    (N.EQ.O)    GO    TO    55 
BRUIT=    AVG/FLOAT(N) 

C  ''"•'    ZI    is    choosen    at    the    intermediate    of    the    five 

C  '•=''    levels    where    the    average    signal    is    stronger. 

\^  ————————————————————————————————————————————— 

C 

IF    (BRUIT. LT.BMIN)    GO    TO    52 

ZI=25.+(DD=M./2./1000.=i=V)    -'■-    FL0AT((I  +  1) +  (IFL-14)) 

LEVEL    =    IFL    +    I    -12 

BMIN    =    BRUIT 

BMAX    =    BIBI 

MAI  =  I 
C 
C 
C 
C 
C 
C 
52  IF    (BRUIT.  GT.BMOU)     GO    TO    55 

BMOU    =    BRUIT 

BMAR(I)    =    BRUIT 
ZL1(D=    25.+(DD=M./2./1000.=:=V)=:=FLOAT((I  +  l)  +  (IFL-14)) 
55  CONTINUE 

C 

C                —    Select    the    value    "ZL"    of    "ZLl"    which    is    below 
C                ■'"''    ZI    but    above    25    meters    which    is    the    first    level. 
C  

c 

DO    62    11=    1,     280 

IF    (ZL1(II).GE.ZI.0R.ZL1(II).LE.25.)    GO    TO    62 
ZL   =    ZLl  (In 
SIG    =    BMAR(II) 
62  CONTINUE 

IF    (BMIN.LE.BBB)    GO    TO    42 
C 
C  '■'•-•'       Compute    the    upper    height    of    the    inversion    layer 

v-/  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  — 

C 

IF    (ZI.EQ.OUT    )    GO    TO    42 
BLIM    =    BMIN    /    2. 
DO    75    II    =    MAI,    NVERT 
AVG    =    0. 
NN       =    0 
DO    70    JJ    =    1^    NZ 

KK=    IFL    +    JJ    +    II    -    2 
IF    (BAC(KK).EQ.OUT    )    GO    TO    70 
NN    =    NN    +    1 
AVG=    AVG    +    BAC(KK) 
70  CONTINUE 

IF    (NN.EQ.O)    GO    TO    75 

70 


c 
c 
c 
c 

r- 

c' 
c 
c 
c 
c 


75 


77 


35 


42 
998 


999 
456 


//GO, 

// 

// 

//GO, 

// 

// 

//GO, 

// 

// 

// 


BRUIT    =    AVG    /    FLOAT    (NN) 
IF    (BRUIT.  LT.BLIM)    GO    TO    77 
SIGU    =    BRUIT 

ZU=25.+(DD='=4./2./1000.=:=V)'!=FLOAT((II  +  l)  +  (IFL-14)) 
CONTINUE 

''""       Estimate    the    thickness    of    the    inversion    layer 

THI    =       ZU    -    ZL 

'■'"''    Write    date,     time,     average  max.     backscatter 

=-:=    signal,     ZI,     SIG,     ZL,     SIGU,     ZU,     THI. 

WRITE(    6_,35)CINT(4J),I=4,6),BMIN,BMAX,ZI,SIG,ZL, 
=i=  SfcU.ZU^THi 

_  WRITE(15  SSHINtr-a,  I) ,  1=4, 6) ,  BMIN,  BMAX,  ZI ,  SIG.  ZL, 

FORMAT    ( IX'  '  r ,'  12 ,  'JUL' ,  IX ,12 ,  'H'I2 , '  T , 

1  F4.0,lX,F4.0,lX,F4.0,'iS4|     ',F4.0,'    ',F4.0,'Mr, 

2  F4.o;'    ',F4.0,'M  ',F4.b,'iI|• 
CONTINU6 

IF    (FLAG.EQ.2)    GO    TO    999 
FLAG    =    2 
GO    TO    111 
WRITE    (6,456) 
FORMAT    (    IX,  58    ('-')) 

STOP 

END 

FTllFOOl    DD    DSN=MSS.S0985.MIZE20,MSVGP=PUB4C, 

DCB=(RECFM=FB_,LRECL=1280,BLKSIZE=6400), 
UNIT=3330V^DISI^=  (OLD, keep! 

FT16F001    DD    DSN=MSS .  ^0985.MIZEi:^,MSVGP=PUB4C, 

DCB=(RECFM=FB_,LRECL=1280,BLKSIZE=6400), 
UNIT=3330V^DISI^=(,OLD,KEEP') 

FT15F001    DD    DSN=MSS  .  ^0985.  ZIBB7,MSVGP=PUB4C, 

DCB=(RECFM=FB_,LRECL=1280,BLKSIZE=6400), 
UNIT=3330V,DISI^=(NEW,CATLG) 
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APPENDIX    B 
TABLE    OF    ZI,    QO,    T* ,    DN/DZ    AND    AM 


[Date    -    Time]     Zi    |    Th    |       CT^     |       Qo       |       T-       jAT    |dN/dZ|ZiM    | 


12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
12Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 


446Z 
454Z 
521Z 
715Z 
723Z 
731Z 
741Z 
749Z 
758Z 

8  7Z 

9  8Z 
917Z 
926Z 
935Z 
943Z 
952Z 

10  2Z 
lOlOZ 

11  2Z 
249Z 
258Z 

8  7Z 
816Z 
834Z 
845Z 

9  3Z 
912Z 
923Z 
932Z 
941Z 
951Z 

10  IZ 
lOlOZ 
1019Z 
1029Z 
1038Z 
1048Z 
1057Z 

11  7Z 
1117Z 
1359Z 
1410Z 
1420Z 
1429Z 
1439Z 
1448Z 
1458Z 

15  8Z 
1517Z 
1526Z 
1536Z 
1545Z 
1555Z 

16  4Z 
1613Z 


493. 

102. 

488. 

88. 

502. 

60. 

474. 

134. 

484. 

125. 

493. 

120. 

502. 

74. 

488. 

144. 

470. 

116. 

470. 

97. 

428. 

70. 

433. 

79. 

424. 

79. 

419. 

65. 

405. 

74. 

386. 

56. 

391. 

60. 

400. 

51. 

396. 

56. 

492. 

132. 

492. 

96. 

449. 

189. 

449. 

195. 

432. 

169. 

406. 

202. 

373. 

109. 

353. 

136. 

339. 

96. 

326. 

96. 

320. 

106. 

313. 

122. 

310. 

116. 

300. 

119. 

356. 

96. 

316. 

132. 

293. 

122. 

326. 

106. 

326. 

93. 

310. 

165. 

234. 

179. 

187. 

89. 

187. 

96. 

197. 

113. 

210. 

126. 

230. 

113. 

234. 

113. 

224. 

99. 

234. 

96. 

217. 

96. 

200. 

122. 

197. 

96. 

210. 

132. 

234. 

139. 

227. 

139. 

217. 

122. 

0011 
0011 
0006 
0016 
0031 
0008 
0017 
0002 
0008 
0005 
0025 
0015 
0012 
0029 
0015 
0011 
0017 
0008 
0029 
0014 
0013 
0024 
0023 
0024 
0016 
0024 
0017 
0032 
0022 
0016 
0015 
0018 
0016 
0019 
0018 
0020 
0020 
0017 
0038 
0028 
0045 
0055 
0036 
0032 
0031 
0032 
0034 
0032 
0033 
0033 
0030 
0036 
0044 
0045 
0049 


0.0019 
0.0026 
0.0047 
0.0014 
0.0012 
0.0016 
0.0014 
0.0014 
0.0020 
0.0015 
0.0038 
0.0027 
0.0028 
0.0025 
0.0033 
0.0053 
0.0048 
0.0048 
0.0067 
0.0016 
0.0011 
0.0013 
0.0013 
0.0012 
0.0010 
0.0013 
0.0011 
0.0017 
0.0016 
0.0014 
0.0021 
0.0025 
0.0032 
0.0040 
0.0048 
0.0043 
0.0051 
0.0051 
0.0032 
0.0032 
0.0022 
0.0217 
0.0258 
0.0207 
0.0141 
0.0141 
0.0170 
0.0118 
0.0119 
0.0121 
0.0119 
0.0115 
0.0097 
0.0095 
0.0081 


0.097 
0.060 
0.081 
•0.019 
0.022 
0.020 
0.027 
0.027 
0.025 
0.050 
0.068 
0.037 
0.038 
0.053 
0.062 
0.111 
0.112 
0.112 
0.299 
0.027 
0.038 
0.022 
0.022 
0.016 
0.012 
0.012 
0.013 
0.018 
0.016 
0.014 
0.019 
0.025 
0.032 
0.033 
0.032 
0.033 
0.028 
0.028 
0.030 
0.041 
0.048 
0.057 
0.063 
0.057 
0.043 
0.043 
0.049 
0.044 
0.042 
0.040 
0.041 
0.040 
0.034 
0.031 
0.031 


0.7 

-7. 

15. 

1.2 

-13. 

13. 

0.4 

-7. 

9. 

5.1 

-38. 

16. 

8.8 

-71. 

11. 

2.6 

-22. 

16. 

4.0 

-55. 

8. 

0.5 

-4. 

22. 

1.9 

-16. 

16. 

0.6 

-7. 

15. 

2.1 

-30. 

9. 

2.3 

-29. 

10. 

1.7 

-22. 

11. 

3.0 

-47. 

7. 

1.3 

-18. 

10. 

0.5 

-9. 

8. 

0.8 

-14. 

9. 

0.4 

-7. 

8. 

0.5 

-9. 

8. 

3.3 

-25. 

17. 

2.1 

-22. 

13. 

6.3 

-33. 

23. 

6.0 

-31. 

25. 

8.8 

-52. 

18. 

7.1 

-35. 

25. 

10.2 

-94. 

7. 

6.6 

-49. 

15. 

8.4 

-88. 

7. 

6.6 

-69. 

8. 

5.5 

-52: 

11. 

3.7 

-31. 

15. 

3.2 

-27. 

15'. 

2.3 

-19. 

16. 

2.8 

-29. 

12. 

2.6 

-20. 

18. 

2.6 

-21. 

17. 

3.4 

-32. 

13. 

2.8 

-30. 

12. 

5.6 

-34. 

20. 

2.6 

-15. 

26. 

3.0 

-34. 

11. 

3.2 

-33. 

12. 

1.9 

-17. 

16. 

1.9 

-15. 

18. 

2.7 

-24. 

15. 

2.8 

-25. 

15. 

2.5 

-26. 

13. 

2.7 

-28. 

12. 

2.8 

-29. 

12. 

2.8 

-23. 

16. 

2.4 

-26. 

13. 

3.2 

-24. 

17. 

4.9 

-35. 

17. 

5.3 

-38. 

17. 

5.7 

-47. 

13. 
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[Date  -  Time  I  Zi  |  Th  |   CT'  |   Qo   |   T'-'   |AT  |dN/dZ|AM 


13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
13Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
14Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 


1624Z 
1633Z 
1642Z 
1653Z 

17  2Z 
1718Z 
1727Z 
1736Z 
1745Z 
1754Z 

18  5Z 
1814Z 
1825Z 
1836Z 
1053Z 
11  OZ 
11  7Z 
1723Z 
1730Z 
1737Z 
1752Z 
1759Z 

18  8Z 
1815Z 
1822Z 
1829Z 
1837Z 
1844Z 
1851Z 
1858Z 

19  6Z 
1913Z 
1920Z 
1927Z 
1942Z 
1949Z 
1957Z 

20  5Z 
228Z 
240Z 
251Z 
316Z 
327Z 
339Z 
352Z 

4  4Z 
415Z 
426Z 
439Z 
451Z 

5  3Z 
514Z 
526Z 
539Z 
550Z 

6  2Z 
615Z 

1036Z 
1049Z 
11  OZ 
1147Z 
1159Z 
1211Z 
1222Z 
1246Z 
1258Z 


234. 

103. 

243. 

109. 

250. 

106. 

260. 

119. 

290. 

152. 

293. 

113. 

296. 

136. 

296. 

96. 

310. 

149. 

323. 

79. 

333. 

129. 

353. 

182. 

363. 

165. 

379. 

83. 

346. 

93. 

366. 

46. 

346. 

129. 

320. 

106. 

320. 

149. 

326. 

89. 

343. 

106. 

343. 

109. 

310. 

103. 

343. 

106. 

333. 

96. 

326. 

79. 

333. 

99. 

330. 

53. 

320. 

103. 

330. 

99. 

320. 

99. 

333. 

106. 

323. 

89. 

330. 

83. 

313. 

103. 

323. 

83. 

330. 

96. 

343. 

93. 

313. 

86. 

310. 

126. 

320. 

142. 

313. 

79. 

310. 

109. 

303. 

126. 

306. 

86. 

303. 

76. 

313. 

139. 

310. 

132. 

310. 

116. 

326. 

106. 

313. 

126. 

303. 

139. 

303. 

129. 

306. 

122. 

313. 

149. 

313. 

113. 

313. 

113. 

286. 

83. 

290. 

79. 

300. 

113. 

300. 

126. 

300. 

86. 

293. 

119. 

293. 

73. 

283. 

126. 

253. 

76. 

0.0045 
0.0016 
0.0056 
0.0012 
0.0022 
0.0017 
0.0009 
0.0010 
0.0013 
0.0021 
0.0012 
0.0008 
0.0014 
0.0011 
0.0009 
0.0004 
0.0005 
0.0018 
0.0009 
0.0011 
0.0011 
0.0018 
0.0018 
0.0011 
0.0033 
0.0042 
0.0025 
0.0022 
0.0015 
0.0018 
0.0024 
0.0018 
0.0013 
0.0038 
0.0013 
0.0026 
0.0027 
0.0014 
0.0040 
0.0055 
0.0033 
0.0027 
0.0049 
0.0025 
0.0030 
0.0033 
0.0040 
0.0032 
0.0036 
0.0037 
0.0025 
0.0042 
0.0029 
0.0027 
0.0036 
0.0014 
0.0017 
0.0030 
0.0030 
0.0044 
0.0048 
0.0047 
0.0042 
0.0049 
0.0015 
0.0033 


0.0069 
0.0058 
0.0040 
0.0023 
0.0015 
0.0047 
0.0068 
0.0077 
0.0036 
0.0043 
0.0045 
0.0042 
0.0033 
0.0024 
0.0012 
0.0012 
0.0011 
0.0012 
0.0015 
0.0015 
0.0014 
0.0016 
0.0016 
0.0024 
0.0017 
0.0017 
0.0017 
0.0023 
0.0025 
0.0025 
0.0023 
0.0014 
0.0011 
0.0011 
0.0015 
0.0015 
0.0014 
0.0011 
0.0016 
0.0010 
0.0015 
0.0020 
0.0027 
0.0038 
0.0043 
0.0038 
0.0022 
0.0034 
0.0032 
0.0023 
0.0021 
0.0028 
0.0032 
0.0020 
0.0017 
0.0020 
0.0018 
0.0013 
0.0011 
0.0010 
0.0012 
0.0018 
0.0043 
0.0036 
0.0048 
0.0048 


0.032 
0.024 
0.021 
0.011 
0.007 
0.017 
0.024 
0.025 
0.015 
0.017 
0.015 
0.014 
0.014 
0.010 
0.014 
0.014 
0.017 
0.027 
0.013 
0.013 
0.015 
0.011 
0.011 
0.014 
0.011 
0.011 
0.011 
0.013 
0.015 
0.015 
0.014 
0.009 
0.006 
0.006 
0.008 
0.008 
0.008 
0.008 
0.005 
0.004 
0.004 
0.006 
0.008 
0.011 
0.013 
0.013 
0.008 
0.011 
0.010 
0.008 
0.008 
0.010 
0.010 
0.007 
0.006 
0.007 
0.006 
0.004 
0.003 
0.003 
0.004 
0.006 
0.014 
0.012 
0.016 
0.017 


5.4 

-53. 

11. 

2.6 

-24. 

15. 

10.9 

-103. 

6. 

4.2 

-35. 

15. 

13.9 

-92. 

10. 

4.3 

-39. 

13. 

1.7 

-13. 

20. 

1.9 

-19. 

13. 

4.0 

-27. 

19. 

5.8 

-73. 

7. 

4.0 

-31. 

16. 

3.0 

-17. 

26. 

5.0 

-30. 

21. 

5.7 

-69. 

7. 

3.3 

-36. 

11. 

1.4 

-30. 

6. 

1.4 

-11. 

19. 

3.0 

-29. 

14. 

3.2 

-22. 

20. 

4.1 

-46. 

10. 

3.3 

-32. 

13. 

8.0 

-74. 

9. 

7.5 

-73. 

9. 

4.0 

-38. 

13. 

14.0 

-146. 

1. 

18.4 

-233. 

-6. 

11.0 

-111. 

5. 

8.2 

-154. 

0. 

4.9 

-48. 

11. 

5.7 

-58. 

10. 

8.0 

-81. 

7. 

9.9 

-93. 

7. 

10.0 

-112. 

4. 

30.0 

-362. 

-17. 

8.2 

-80. 

8. 

16.2 

-195. 

-3. 

15.7 

-163. 

-1. 

8.2 

-89. 

6. 

35.6 

-413. 

-2'? 

68.8 

-546. 

-49! 

34.3 

-242. 

-12. 

21.2 

-268. 

-9. 

26.9 

-247. 

-10. 

10.  1 

-80. 

10. 

10.7 

-124. 

3. 

11.4 

-149. 

1. 

22.4 

-161. 

-1. 

13.2 

-100. 

8. 

16.5 

-142. 

2. 

23.0 

-217. 

-6. 

14.6 

-116. 

5. 

18.5 

-133. 

3. 

12.7 

-99. 

8. 

18.1 

-148. 

1. 

27.4 

-184. 

-4. 

9.6 

-85. 

8. 

12.4 

-110. 

5. 

29.9 

-361. 

-17. 

39.9 

-505. 

-28. 

57.0 

-504. 

-39. 

52.8 

-419. 

-33. 

34.8 

-405. 

-21. 

13.1 

-110. 

6. 

17.7 

-243. 

-6. 

4.2 

-33. 

16. 

7.7 

-102. 

4. 
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Date    -    Timel    Zi    |    Th    |       CT^    |       Qo       |       T-       |  AT    |dN/dZ|AM    | 


15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
ISJul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 
15Jul 


1310Z 
1323Z 
1334Z 
1345Z 
1359Z 
1410Z 
1422Z 
1433Z 
1444Z 
1457Z 

15  9Z 
1521Z 
1532Z 
1545Z 
1557Z 

16  8Z 
1620Z 
1737Z 
1748Z 
18  OZ 
1811Z 
1823Z 
1836Z 
1847Z 
1858Z 
1912Z 
2010Z 
2022Z 
2034Z 
2047Z 
2059Z 
2111Z 
2122Z 
2134Z 
2147Z 


260. 

106. 

273. 

113. 

273. 

93. 

277. 

93. 

283. 

96. 

286. 

86. 

283. 

96. 

293. 

116. 

303. 

83. 

306. 

129. 

316. 

76. 

320. 

73. 

316. 

139. 

306. 

96. 

296. 

89. 

293. 

132. 

263. 

106. 

250. 

93. 

253. 

89. 

253. 

93. 

260. 

106. 

277. 

116. 

277. 

116. 

277. 

113. 

263. 

73. 

260. 

96. 

243. 

79. 

257. 

89. 

263. 

89. 

263. 

96. 

267. 

99. 

270. 

99. 

263. 

106. 

270. 

109. 

277. 

106. 

0028 
0027 
0026 
0049 
0033 
0031 
0032 
0025 
0038 
0029 
0018 
0025 
0013 
0023 
0040 
0028 
0022 
0025 
0033 
0048 
0020 
0015 
0032 
0018 
0044 
0041 
0038 
0063 
0033 
0055 
0023 
0032 
0028 
0042 
0020 


0.0051 
0.0053 
0.0066 
0.0039 
0.0050 
0.0033 
0.0031 
0.0020 
0.0014 
0.0030 
0.0035 
0.0035 
0.0039 
0.0040 
0.0053 
0.0038 
0.0037 
0.0060 
0.0069 
0.0065 
0.0067 
0.0066 
0.0056 
0.0068 
0.0065 
0.0052 
0.0051 
0.0025 
0.0033 
0.0031 
0.0030 
0.0029 
0.0020 
0.0013 
0.0016 


0.014 
0.016 
0.017 
0.011 
0.014 
0.011 
0.008 
0.005 
0.004 
0.009 
0.010 
0.013 
0.014 
0.016 
0.016 
0.014 
0.014 
0.018 
0.023 
0.021 
0.022 
0.024 
0.019 
0.021 
0.021 
0.016 
0.015 
0.009 
0.011 
0.010 
0.010 
0.009 
0.006 
0.004 
0.005 


8.2 

-77. 

8. 

7.1 

-63. 

11. 

6.5 

-70. 

8. 

18.7 

-201. 

-4. 

9.8 

-102. 

5. 

11.7 

-136. 

2. 

16:7 

-174. 

-2. 

21.6 

-186. 

-3. 

37.8 

-455. 

-25. 

14.4 

-112. 

6. 

7.9 

-105. 

4. 

8.8 

-121. 

3. 

4.3 

-31. 

18. 

6.4 

-67. 

9. 

10.7 

-121. 

3. 

9.0 

-68. 

12. 

6.3 

-59. 

10. 

5.5 

-59. 

9. 

5.7 

-64. 

8. 

8.9 

-96. 

6. 

3.7 

-35. 

13. 

2.7 

-23. 

16. 

7.1 

-62. 

11. 

3.5 

-31. 

14. 

8.7 

-119. 

3. 

10.5 

-110. 

5. 

9.7 

-123. 

3. 

28.6 

-321. 

-15. 

12.3 

-138. 

2. 

23.4 

-244. 

-8. 

9.5 

-96. 

6. 

15.6 

-158. 

0. 

17.8 

-168. 

-1. 

45.1 

-414. 

-28. 

18.1 

-171. 

-1. 

74 


LIST    OF   REFERENCES 


Businger,  J.  A.  (1973):  Turbulent  transfert  in  the  atmospheric 
surface  layer.  Workshop  on  Micrometeorology .  D.A.  Haugen  (ed.), 
American    Meteorol.     Soc,     Boston,     pp    1-65. 

Davidson,  K.L,  T.M.  Houlihan,  C.W.  Fairal  and  G.E.  Schacher 
(1978):  Observation  of  the  temperature  structure  function  param- 
eter, CT^ ,  over  the  ocean.  Boundary-Layer  Meteor.,  15,  pp 
507-523. 

Dorman,  C.E.  (1985):  Evidence  of  Kelvin  waves  in  California's 
marine  layer  and  related  eddy  generation.  .  Monthly  weather  review, 
113,    pp. 827-839. 

Fairal,  C.W.,  R.  Markson,  G.E.  Schacher  and  K.L.  Davidson 
(1980):  An  aircraft  study  of  turbulence  dissipation  rate  and 
temperature  structure  function  in  the  unstable  atmospheric  boundary 
layer.        Boundary-Layer    Meteorol.,     19,     pp    453-469. 

Fleagle,  R.G.  and  J.  A.  Businger  (1980):  An  Introduction  to 
Atmospheric  Physics.  second  editon.  Academic  Press.  N.Y.,  432 
pp. 

Friehe,  C.A.  (1977):  Estimation  of  the  refractive-index  temperatut^e 
structure    parameter    over    the    ocean.       App.     Opt.,     16,     pp.     334-340. 


75 


Hitney,  H.V.,  R.A.  Paulus,  C.P.  Hattan,  K.D.  Anderson  and 
G.E.  Linden  (1981):  IREPS  Revision  2  User's  Manual.  Pub. 
Naval    Ocean    Systems    Center.    San    Diego,     CA.,     134    pp. 

Johannessen,  O.M.  and  D.A.  Horn  (1984):  MIZEX  84-  Summer 
experiment   preliminary   report.       CRREL    Special    report    84-29. 

Meeks,  M.L.  (1984):  Radar  Propagation  at  Lo\\>  Altitude.  Artech 
house,    inc.    Dedham,    MA.,     105    pp. 

Mohn,  M.  (1985):  The  structure  of  the  atmospheric  boundary  layer 
in  the  marginal  ice  zone  as  derived  from  sodar.  (NPS  thesis  in 
preparation) . 

Monin,  A.S.  and  A.M.  Obukhov  (1954):  Basic  Zaws  of  turbulent 
mixing  in  the  ground  layer  of  the  atmosphere.  Akad.  Nauk  SSSR 
Geofiz.     Inst.    Tr.,     151,    pp.     163-187. 

Parry,  D.H.  and  Melvin  J.  Sanders,  Jr.  (1972):  The  design  and 
operation  of  an  acoustic  radar.  IEEE.  Trans,  on  Geosci.  Elect., 
10,    pp.    58-64. 

Parry,  H.D.,  Melvin  J.  Sanders.  Jr  and  H.P.  Jensen,  (1975): 
Operational  applications  of  a  pure  acoustic  sounding  system.  J. 
Appl. Meteor.,     14,     pp.     67-77. 

Wyngaard,  J.C.,  (1973):  On  surface  layer  turbulence.  Chapter  3, 
Workshop  on  micrometeorology,  D.  Haugen,  Editor,  American 
Meteorological    Society,     pp.     101-149. 


76 


Wyng-aard,  J.C.,  Y.  Izumi  and  S.A.  Collins  (1971):  Behavior  of 
the  refractive  index  structure  parameter  near  the  ground.  J.  Opt. 
Soc.     Am.,     61,     pp    1646-1650. 

Wyngaard,  J.C.  and  M.A.  Lemone,  (1980):  Behavior  of  the 
refractive -index  structure  parameter  in  the  entraining  convective 
boundary    layer.       J.    Atm.    Sci.,    37,    pp.    1573-1585. 


77 


INITIAL   DISTRIBUTION   LIST 


No.       Copies 


1.  Defense    Technical    Information    Center  2 
Cameron    Station 

Alexandria,    Virginia    22304-6145. 

2.  Library^     Code    0142  2 
Naval    Postgraduate    School 

Monterey,     California    93943. 

3.  Professor    C.N.K.    Mooers,     Chairman  1 
Department    of    Oceanography,     code    68Mr 

Naval    Postgraduate    Sciiool 
Monterey,     C^alifornia    93943. 

4.  Professor    R.J.     Renard,     Chairman  1 
Department    of    Meteorology,     Code    63Rd 

Naval    Postgraduate    School 
Monterey,    California    93943. 


5.  Professor    W.J.     Shaw 

Department    of    Meteorology, 
Naval    Postgraduate    School 
Monterey,     California    93943. 


code    63 


Professor    R.W.     Garwood 

Departement    of    Oceanography,     code    68 

Naval    Postgraduate    School 

Monterey,    California    93943. 


7.  Professor    K.L.     Davidson 

Department    of    Meteorology,     code    63 
'    ""     tgraduate    School 


Naval    Pos 

Monterey,    California    93943. 


8.  P.     Guest 

Department    of    Meteorology,     code    63Gs 
Naval    Postgraduate    School 
Monterey,    California    93943. 


Doctor    Leonard    Johnson 
Office    of    Naval    Research 
Arctic    science 

800    N.     Quincy    St.,     BCT    No    1, 
Arlington,     VA    22217-5000. 


10.        Doctor    Paul    TwitcheU 
Code    Air-330 

Naval    Air    Systems    Command, 
Washington    D.C.     20360. 


78 


11.  Docteur    J.C.    Gascard 
Laboratoire    d'Oceanographie    Physique 
Museum    d'Histoire    Naturelle 

43.     rue    Cuvier 

75z31  Paris  cedex  05,  France, 

12.  Docteur    A.    Weill 

Centre    de    Recherche    en    Physique    de    I'Environnement 

CNET 

38-40,     rue    du    General    Leclerc 

92131    Issy-les   -Moulineaux,     France. 


13.        Sheryl    Fellbaum 
Tamar    Neta 
Pat    Boyle 

Department    of    Meteorology, 
Naval    Postgraduate    School 
Monterey,     California    93943 


code    63Ds 


14.  Monsieur    I'Amiral 

Chef    d'Etat-Major    de    la    Marine 

2,     rue    Royale 

75200    Paris    Naval,     France. 

15.  Monsieur    le    Capitaine    de    Vaisseau 

Attache    Naval    aupres    de    I'Ambassade    de    France 
4101,     Reservoir    road    NW 
Washington    D.C.     20007. 

16.  Monsieur    le    Capitaine    de    Vaisseau 
Commandant    le    CEPMAN 

B.P.     38,     83800    Toulon    Naval,     France. 

17.  Capitaine    de    Vaisseau    (e.r)     Beydon 
Meteorologie    Nationale    (Meteo/Maritime) 
92106    Boulogne    Billancourt    Cedex,     France. 

18.  Monsieur    le    Directeur    de    I'EPSHOM 
13,     rue    du    Chatellier 

B.^.     426 

29275    Brest    Cedex,     France. 

19.  Capitaine    de    Corvette    J.F.     Soubrier 
Etat-Major    de    la    Marine    a    Lorient 
B.P.     467 

56324    Lorient,     France. 

20.  Lieutenant    de    Vaisseau    E.     de    Rouge 

Cellule    Operationelle    d'Environnement    de    la    Marine 
2,    avenue    Rapp 
75007    Paris,     France. 

21.  Meteorologie    Nationale 
SMM    /    Documentation 
2,    avenue    Rapp 
75007    Paris,     France. 


79 


i^q  -^/? 


21i^l5S 


LL 


8 


De   Rouge 

An  investigation  of 
the  atmospheric  boun- 
dary layer  over  the 
Artie  Ocean  using 
SODAR. 

JLC    tG  i  i  6  ti 


15  APR   87 

21  OCT    ^'7 

21  OCT  a? 


I 


333II0 

S  1  05  9  1 


Thesis 

D4424 

c.l 


2:a:5s 


De  Rouge 

An  investigation  of 
the  atmospheric  boun- 
dary layer  over  the 
Artie  Ocean  using 
SODAR. 


